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A general field-theoretical description of many-fermion systems, both with and without quenched 
disorder, is developed. Starting from the basic Grassmannian action for interacting fermions, we 
first bosonize the theory by introducing composite matrix variables that correspond to two-fermion 
excitations and integrating out the fermion degrees of freedom. The saddle point solution of the 
resulting matrix field theory reproduces a disordered Hartree-Fock approximation, and an expansion 
to Gaussian order about the saddle point corresponds to a disordered RPA-like theory. In the 
clean limit they reduce to the ordinary Hartree-Fock and random-phase approximations. We first 
concentrate on disordered systems, and perform a symmetry analysis that allows for a systematic 
separation of the massless modes from the massive ones. By treating the massive modes in a simple 
approximation, one obtains a technically satisfactory derivation of the generalized nonlinear sigma- 
model that has been used in the theory of metal-insulator transitions. The theory also allows for the 
treatment of other phase transitions in the disordered Fermi liquid. We further use renormalization 
group techniques to establish the existence of a disordered Fermi-liquid fixed point, and show that 
it is stable for all dimensions d > 2. The so-called weak-localization effects can be understood as 
corrections to scaling near this fixed point. The general theory also allows for studying the clean 
limit. For that case we develop a loop expansion that corresponds to an expansion in powers of 
the screened Coulomb interaction, and that represents a systematic improvement over RPA. We 
corroborate the existence of a Fermi-liquid fixed point that is stable for all dimensions d > 1, in 
agreement with other recent field-theoretical treatments of clean Fermi liquids. 
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I. INTRODUCTION 

The many-fermion problem has a long history, due 
to its importance with respect to electrons in condensed 
matter systems.ETu In recent years there has been a re- 
newed interest in fundamental aspects of this problem, 
especially for electrons at low or zero temperature, both 
with and without quenched disorder. The important 
physical systems for which these studies are relevant in- 
clude, superconductors, both high-T c and conventional 
ones, doped semiconductors in both their metallic and 
insulating phases, amorphous alloys, and Quantum Hall 
systems. 

Historically, there have been two important techniques 
to tackle the many-fermion problem: Landau's phe- 
nomenological Fermi-liquid theoryQ and the microscopic 
many-body-, perturbation theory or Feynman diagram 
approach,0 : which starts from a canonically quantized 
Hamiltonian of the problem. While both techniques are 
very useful, and have yielded many important results, 
either one of them has serious limitations. Fermi-liquid 
theory has been extended to include disorder ,Eli3 but often 
a microscopic approach is preferred. For the microscopic 
many-body perturbation theory, the inclusion of disorder 
is extremely awkward, mostly because of a large num- 
ber of diagrammatic contributions that are individually 
finite, but that ultimately cancel each other. More re- 
cent approaches have used functional or field-theoretical 
methods,0 which allow for the application of renormaliza- 
tion group (RG) techniques. A significant recent advance 
in this area has been Shankar's RG technique^ for clean 
fcrmion systems that is applied directly to the Grassman- 
nian field theory for fermions. This approach, which is 
still in its infancy, is very promising. For instance, it has 
already led to a derivation of -Fermi-liquid theory from 
a microscopic starting point. Em It also has been used to 
provide RG derivations of the Cooper instability prob- 
lem, and of RPA-like screening.!! Other applications, for 
instance a discussion of corrections to scaling near the 
RG fixed point (FP) that describes the clean Fermi liq- 
uid, should be possible. However, for reasons that will 
be discussed in detail in the present paper, it is not easy 
to include the effects of quenched disorder into this ap- 
proach either. 

An alternative approach for clean electronic systems 
has been pursued by Houghton, Marston, and others,E2l 
and by Frohlich and Gotschmann.EJ These authors have 
generalized the bosonization techniques that have been 
used successfully in d = 1 to higher dimensions. In Ref. 
|TT] the bosonization was also combined with RG tech- 
niques. These theories have been used, for example, 
to rederive certain nonanalyticities that occur in Fermi- 
liquid theory, and to study the possibility of marginal 



Fermi liquids in dimensions d > 2. Finally, clean 
electronic systems have been bosonized by means of a 
Hubbard-Stratonovich transformation and the use of the 
classical HubbardpStratonovich field as the fundamental 
field of the theory.lia Within this latter approach, disorder 
has been included in the single-particle Green function at 
the level of the lowest order Born approximation.^ More 
sophisticated effects of quenched disorder have so far not 
been included into any of these approaches. 

For systems of disordered electrons, a completely dif- 
ferent theory has been in use for some time, that de- 
scribes composite fermionic, i.e. effectively bosonic, de- 
grees of freedom. For noninteracting electrons, this ef- 
fectivep&eld theory takes the form of a nonlinear sigma- 
model,li-3 which was generalised to the case of interacting 
electrons by Finkel'stein.ll3ll3 It is custom tailored for the 
description of the metal-insulator transition near d = 2, 
and it does not allow for the clean limit to be taken. No 
technically satisfactory derivation of the interacting the- 
ory has ever been given, but rather its structure has been 
guessed, based on more rigorous derivations of Wegner's 
effective model for the metal-iswilator transition in non- 
interacting electronic systems .113113 The key idea underly- 
ing these effective field theories is to keep explicitly only 
those degrees of freedom that are likely to be relevant for 
the problem under consideration, and to integrate out all 
others in some simple approximation. 

It is clearly desirable to develop a more flexible field- 
theoretical approach for the many-fermion problem. 
Some of the desired features of such a theory are as fol- 
lows. (1) Both clean and disordered systems should be 
describable within a single framework. (2) The theory 
should allow for a RG description of both clean and disor- 
dered Fermi liquids, of the metal-insulator transition, as 
well as of other quantum phase transitions, like e.g. the 
recently investigated magnetic ones.lia~Eil (3) It should 
offer flexibility with respect to how many, and which, de- 
grees of freedom one keeps explicitly, depending on the 
problem under consideration. (4) It should provide a sat- 
isfactory derivation of the nonlinear sigma-model used 
to describe the metal-insulator transition in interacting 
disordered electronic systems. (5) It should allow for ex- 
plicit calculations of physically relevant observables, such 
as thermodynamic and transport properties. While all of 
the existing theories fulfill some of these criteria, none of 
them meet all of them simultaneously. 

This is the first of two papers where we develop such a 
theory. In the present paper we discuss our general phys- 
ical ideas, and mainly focus on the disordered interacting 
fermion problem, although we will discuss some aspects 
of the clean limit. In a future, second paper we plan 
to use the same approach to thoroughly discuss clean 
interacting fermion systems, as well as the connection 
between our approach and others. The outline of this 
paper is as follows. In Sec. [j] we introduce our start- 
ing point, a microscopic model for an electron fluid, in 
general in the presence of quenched disorder. We then 
show how to transform this Grassmannian field theory 
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into one for classical or bosonic fields, and discuss the 
physical motivation for this transformation. In this sec- 
tion we also construct a saddle-point solution for the 
resulting composite-fermion field theory, and expand to 
Gaussian order about the saddle point. The result is 
equivalent to what one obtains within many-body di- 
agrammatic theory from the random phase approxima- 
tion (RPA), modified by disorder. In Sec. Ill we perform 
a symmetry analysis of our field theory. The results, to- 
gether with some conclusions drawn from Sec. |n[ suggest 
slightly separate ways to proceed from here for clean and 
disordered fermions, respectively. Focusing on the latter, 
we use the symmetry analysis to identify and classify all 
of the slow or soft modes of the system. This is crucial 
for a RG description of the problem. Using these results, 
we then derive an effective field theory for the slow modes 
in a disordered fermion system. We identify a FP that 
describes a disordered Fermi liquid, and show that it is 
stable for d > 2. We then show that the so-called weak- 
localization effects can be interpreted as corrections to 
scaling near this FP. We derive Finkel'stein's generaliza- 
tion of Wegner's nonlinear sigma-model that has been 
used to describe metal-insulator transitions near d = 2, 
and discuss the critical FP that corresponds to this quan- 
tum phase transition. The section is concluded by a 
brief discussion of a magnetic phase transition that is 
described by the theory, and of the relation between the 
theory presented here and earlier theories of magnetic 
transitions in disordered systems. In Sec. IV we consider 



the clean limit of our theory and discuss the clean Fermi- 
liquid FP. There we show that certain nonanalyticities 
that appear in Fermi-liquid theory can be interpreted as 
corrections to scaling near this FP, and we point out far- 
reaching analogies between clean and disordered Fermi 
systems. In Sec. ^ we present our conclusions. Appendix 
A spells out a technical point that one encounters in Sec. 
, and Appendix |b] contains a pedagogical discussion of 
an O(N) symmetric <j> theory. 



II. MATRIX FIELD THEORY 

A. Grassmannian field theory 

Let us start with a field-theoretical description of a 
system of interacting- disordered fermions in terms of 
Grassmann variables.Q The partition function of the sys- 
tem is 



Z= / D[^] e 



sbl>,i>] 



(2.1) 



Here the functional integration is with respect to Grass- 
mann valued fields, tp and ip, and the action, S 1 , is given 
by 

S = - ^Pcr(x) d T i)) a (x) + S + S dls + S int ■ 



We use a (d + l)-vector notation, with x — (x, r), and 
J dx = J v dx J dr. x denotes position, r imaginary 
time, V is the system volume, f3 = 1/T is the inverse 
temperature, a is the spin label, and we use units such 
that H = ks = 1- 5*o describes free fermions with chemi- 
cal potential jj,, 



S = J dx^ Vv(a 



2m 



fj, lp a (x) 



(2.2b) 



with m the fermion mass. The Laplacian will be denoted 
by V 2 throughout. We will mostly be concerned with 
systems at T = 0, where /i = ep, with ep = k F /2m 
the Fermi energy, and hp the Fermi momentum. Sdis 
describes a static random potential, it(x), coupling to 
the fermionic number density, 



S, 



dis 



J dx^2 u{x)il} rT {x)ijj a {x) 



(2.2c) 



and Si„t describes a spin-independent two-particle inter- 
action, 



Si n t = ~ J dxx dx 2 ^2 u(xi-x 2 ) 



xil) a - 1 (x 1 )ip a2 (x 2 )ip l T 2 (x2)tpa- 1 (xi) . (2. 2d) 

The interaction potential i>(x) will be specified below . 

For simplicity, the random potential in Eq. ( 2.2c ) is 
taken to be Gaussian distributed with a second moment 
that is given by a function U (x) with the dimension of 
an energy density, 



{u(x)u(y)} dis = 



1 



nNp 



^(x-y) 



(2.3) 



where {. . .} dis denotes the disorder average. Np is a 
normalization factor with the dimension of a density of 
states. One could use the free electron density of states at 
the Fermi level, but it will turn out to be more convenient 
to include some trivial disorder and interaction renormal- 
izations. We will explicitly define Np in Eq. ( 2.46| ) below. 
The full physical density of states at the Fermi level will 
also be encountered later, and will be denoted by N(ep). 
More general disorder potentials can be considered, but 
as long as they couple only to the electron number den- 
sity, the diffe rence between the more general potentials 
and Eq. ( [2.3|) can be shown to be RG irrelevant at the 
disordered Fermi-liquid FP. For other purposes, e.g. a 
description of the Anderson-Mott metal-insulator tran- 
sition in the system, the difference is likewise expected 
to be RG irrelevant. 

Since the system contains quenched disorder, it is nec- 
essary to average the free energy or Jn.Z. This is accom- 
plished by means of the replica trickjE 2 ] which is based on 
the identity 



(2.2a) 



InZ 



lim (Z r ' 



l)/n 



(2.4) 
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Introducing n identical replicas of the system (with n an 
integer), labeled by the index a, and carrying out the 
disorder average, we obtain 



p n 

Z = {Z n } dis = J ]\D exp[£] 



(2.5) 



We again separate S into free, disordered and interaction 
parts as follows, 



n 



(2.6a) 



It is useful to go to a Fourier representation with wave 
vectors k and fermionic Matsubara frequencies u> n — 
1-KT(n + 1/2), and a (d+ l)-vector notation, k = (k, u> n ). 
Then 

So = E <W K - k V2m + M ] , (2.6b) 

and, 



2%N F 



i+k :i ,k2+k4 



tnz EE EE 

/3=1 {ki} n,m o,ct' 

x[/(k 1 - k 2 ) ^(k 1 )C.(k 2 )^(k 3 ) 

xV^ iCT ,(k 4 ) , (2.6c) 



T 
2" 



E E S ki+k 2 ,k 3 +k4 u ( k 2 -k 3 ) 
0-1,0-2 {ki} 

xC.^iXffeXffeli"^) • (2.6d) 
Here we have defined 

i> na (k) = V^ CT (A) = V^A 7 / e - 4 ( kx — Vv (a:) , 



(2.7a) 



Vw(k) = Vv(fc) = / dx e l( - 



(kx- 



(2.7b) 

Ultimately, we will be interested in long- wavelength, 
low-frequency processes. For clean electrons, this in gen- 
eral means that only the scattering of electrons and holes 
close to the Fermi surface is important. In a clean Fermi 
liquid this is true because the lifetimes of the single- 
particle excitations become infinitely long as the Fermi 
surface is approached. For the disordered electron prob- 
lem the situation is different. All single-particle momen- 
tum and frequency excitations are damped by the single- 
particle relaxation rate, l/r re i, which in general is not 
small. For this case, the only slow or soft excitations 
are two-particle excitations. This implies that the soft 
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FIG. 1. Typical small-angle (1), large-angle (2), and 
2&F-scattering processes (3) near the Fermi surface in d = 2. 

excitations in the two cases are fundamentally different, 
and that clean systems have in general more soft modes 
than disordered ones. This observation will be of great 
importance below, and will suggest somewhat different 
effective field theories for dealing with clean and disor- 
dered systems, respectively. In fact, this^difference is the 
reason why the approaches of Shankara and Houghton 
and MarstonE2l cannot be simply extended to the dirty 
case. 

To examine this important point in more detail, we 
first consider the clean case. It is customary and con- 
venient to divide the possible scattering processes into 
three classes: (1) small-angle scattering, (2) large-angle 
scattering, and (3) 2fcp-scattering. These classes are also 
referred to as the particle-hole channel for classes (1) and 

(2) , and the particle-particle or Cooper channel for class 

(3) , respectively. The corresponding scattering processes 
are schematically depicted in Fig. [y. We next make ex- 
plicit the phase space decomposition t hat is inherent to 
this classification by writing S^ t , Eq. ( 2.6d ), as 



,a(2) 



ia(3) 



(2- 



where, 



5«(1) 

J int 



~2 EEE v (q) c,w«>+?) 

0"i<T2 k,p q 

x<,G#°(fc + «) , (2.8b) 



S 



a (2) 



EEEVp- 

(T\(T2 k,p Q 



k) CWC(|) + 9) 
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xC 2 (fc + <z)C» 



(2.8c) 



s?J 3) = -^ E EE' w ( k +p)<(- fc X(fc+9) 



x< 2 (-p + ?)C» 



(2.8d) 



Here the prime on the g-summation indicates that only 
momenta up to some momentum cutoff, A, are integrated 
over. This restriction is necessary to avoid double coun t - 
ing, since each of the three expressions, Eqs. ( 2.8b ) - 
(2.8d), represent all of S^ t if all wave vectors are summed 
over. In general the long- wavelength physics we are in- 
terested i n wi ll not depend on A. Finally, one often writes 
the Eqs. (2.8) in terms of singlet (s) and triplet (t) con- 
tributions. We introduce spinors 



and their Fourier transforms 



(2.9a) 



(2.9b) 



as well as their adjoints i V> a (fc), and a scalar product in 
spinor space, (ip,ip) = tp ■ ip, where the dot denotes the 
matrix product. Then we can write the interaction term 
as 



qa _ qa(s) q a(t) qa(3) 
J int — °int ' °int ' °int i 



(2.10a) 



with, 



S 



(2.10b) 



S 



a (t) 



k,p q 



A:,p <? 2—1 

x(V a (p + g),«iV a (p)) 



(2.10c) 



Here Sj = icTj, with <Ti,2,3 the Pauli matrices, and sq = <7q 
is the 2x2 identity matrix. We have also defined the 
singlet (s) and triplet (t) interaction amplitudes 



r£(<?) = ^(p-k) , 



(2.11a) 



Th e effec tive i nteraction potentials that appear in Eqs. 
( |2.10| ) and ( 2.11 ) are all given by the basic interaction po- 
tential v, taken at different momenta. S^j: 1 ^, and S?j: 
as well, contains the direct scattering contribution, or 
v(q), with q the restricted (|q| < A << kp) momentum. 
If v is chosen to be a bare Coulomb interaction, then 
the direct scattering contribution leads to singularities 
in perturbation theory that indicate the need for infinite 
resummations to incorporate screening. For simplicity, 
we assume that this procedure has already been carried 
out, and take v to be a statically screened Coulomb inter- 
action. For sm all |q|, it is then sufficient to replace i>(q) 
in Eq. ( 2.11 b| ) by a number. The moduli of the other 
momenta, k and p in Eqs. (2.10), are equal to kp for the 
most important scattering processes since, as mentioned 
above, the soft or slow excitations in the clean Fermi 
system involve particles and holes close to the Fermi sur- 
face. The angular dependences of these coupling con- 
stants is usually taken into account by expanding them 
in Legendre polynomials on the Fermi surface a For fu- 
ture reference we note that if only the zero angula r mo - 
mentum channel, I = 0, is retained, then Eqs. ( |2.10| ) 

and (|I|) are valid with rjg,(0), rg(0), and rg(0) 
replaced by numbers, T^ s \ and T^ c \ which are the 
Fermi surface averages of the respective interaction am- 
plitudes. Shankar,B in his RG approach, has explicitly 
derived these results, and has shown that all corrections 
to these approximations are RG irrelevant with respect 
to a RG fixed point that describes a Fermi liquid. 

For disordered Fermi systems the situation is different, 
since single-particle excitations are massive modes with 
a mass that is proportional to the inverse elastic mean- 
free time. In Section III we will give a detailed discussion 
of the soft or slow modes in a disordered Fermi system. 
The conclusion will be that, of the modes that appear 
in Si n t, the dominant soft modes are those that involve 
fluctuations of either the particle number density, n n , or 
the spin density, n s , in the particle-hole channel, or dee-, 
sity fluctuations, n c , in the particle-particle channel.cij 
I n ter ms of our spinors, the latter are given as (cf. Eqs. 



(2.12a) 



<(<?) = v/T/yE^W > s ^ a ( k + l)) . (2.12b) 



and 



rg(<?) = v(q)-^(p-k) 



(2.11b) 



In addition we define the Cooper channel or 2/cf- 
scattering amplitude 



r£>(g) = «(k + P ) 



(2.11c) 



(2.12c) 



In Eqs. ( 2.10|) this implies that we should project onto 
these density fluctuations. This projection is performed 
in Appendix The net result is that the interaction 
amplitudes can be replaced by constants, namely their 
angular averages with a distribution of momenta near 
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the Fermi surface. We thus obtain the interaction part 
of the action in the form, 

^ (S) = -\ r (s) E E ' (^ a ( fc )' s or(k + q )) 

k,p q 

x(^> + g),s V») , (2.13a) 

= -\ rw E E ' £dr (* w-(* + <?)) 

x(V°(p + g),SiV a (p)) • (2.13b) 



5V 



<(3) _ 



<Ti/<T 2 fc,p g 
xC 2 (P + 9)i"("P) 



C x (fc)C 2 (-fc + 9) 



(2.13c) 



The disordered Fermi system is thus described by only 
three interaction p aram eters, r( s ), rw, and r^ c ', which 
are given by Eq. (A4b) and analogous expressions. As 
long as we are only interested in physics that is dom- 
inated by soft modes, this restriction is exact. Clean 
Fermi systems, on the other hand, in general require an 
infinite number of Fer mi-li quid parameters, and the re- 
striction to the Eqs. ( 2.13 ), with the three interaction 
constants related to the Fermi-liquid parameters in the 
£ = channel, constitutes an approximation. Physically 
this difference is due to the fact that there are more soft 
modes in the clean case than in the dirty case. 



B. Composite variables: Matrix field theory 

1. The action in terms of Q -matrices 

As noted above, the slow modes for the disordered 
fermion problem are fluctuations of products of fermion 
operators. It is therefore convenient to transform to a 
field theory in terms of these composite variables. For 
both technical and physical reasons, which will become 
clear in Sec. HI, it is convenient to go to a bispinor and, 
eventually, to .a_spin-quaternion representation. We de- 
fine a bispinorn.3 



±f <(x) 
V2 V s 2^(x) 



1 

71 



<(x) 

V-C T «/ 



and an adjoint bispinorEl 



(T7+)S(x) = i<(x) = -j= 



/-V£ T (x)\ 

Ci(x) 
V-<(x)/ 



(2.14a) 



(2.14b) 



with c the charge-conjugation matrix 



s 2 
s 2 



IT\ <g> S 2 



(2.14c) 



The four degrees of freedom represented by the bispinor 
are the particle-hole or number density degrees of free- 
dom, and the two spin degrees of freedom. We have 
also defined a basis in spin-quaternion space as r r ® s, 
(r, z = 0,1, 2, 3), with t = Sq the 2x2 identity matrix, 
and Tj — —Sj — —i(Jj (j — 1,2,3), with Gj the Pauli 
matrices. In this basis, the channels r = 0, 3 and r = 1, 2 
describe the particle hole and particle-particle degrees of 
freedom, and the channels i = and i = 1,2,3 describe 
the spin-singlet and spin-triplet, respectively. 

It is convenient to define a scalar product in bispinor 
space as 



( , )»°.it)=IijC , ')! 1 » 



(2.15) 



where J rj denotes the elements of 77, and ^ 77 denotes the 
elements of 77 + . The adjoint of any operator A in bispinor 
space with respect to this scalar product is given by 
c T A T c. In terms of the bispinors, the terms on the right- 
hand side of Eq. (|2.6a| ) can be written 



So - -*E ( ? > Q ( fc )' ^« - k2 / 2m + d ^ a ( fc )) 



(2.16a) 



S, 



dis 



and 



2ttN, 



EEE 5 ki+k 3 ,k 2 +k 4 U (ki - k 2 ) 



a, {3 n,m {ki} 

x«(k 1 ),<(k 2 ))(^(k 3 ),^(k 4 )) , (2.16b) 



c c ( s ) 1 c W 1 c I 

>->mt — t~>i nt 1 J int J v, 



(2.17a) 



5(«) 



EEE E 

a k,p q r=0,3 



i-iy 



x(r ? Q (fc),(r I .® So )r 7 a (fc + g )) 

x(r 1 a (p + q)ATr®s )r 1 a (p)) , (2.17b) 



S 



■i1~EEE'EM)' 



3 

E 

i=l 



a k,p q r=0,3 

x(ti a (k),(T r ®8i)rj*(k + q)) 

x((r) a (P + q)ATr®s l )r ] a (p)) , (2.17c) 



5(c) _ 

°int ~ 



J>p(c) 



EEE' 



, E E 

a k,p q ni ,n 2 ,m r=l,2 



x «(-k), (r r <g> s ) r,« Ili+m (-k + q)) 
x (r£„ 3 (-p), (r r ® s ) < 2+m (-p - q)) 



(2.17d) 
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The same phase space decomposition used to break up 
Si n t can also be used to rewrite Sdis as 



5 _ 5(1) , 5(2) 
Jdis — O dis ~r O dig , 



(2.18a) 



with, 



5i2 = 



(2.18b) 



* S 27TiV F Tl 

q,/3 n,m k,p q 
X(??m(p),»7m(p-q)) » 



5(2) _ 



1 



EEEE'te( k ).^(p)) 



a,/3 n,rn k,p q 

x(^(p + q),^(k + q)) 



(2.18c) 



Here we have introduced two different scattering times, 
ti and T re i, that will in general renormalize differently. 
The reason for their appearance is related to the poten- 
tial correlation function U(ki — k^) in Eq. (2.161), which 



leads to a U(q) in sj^, and a U(k — p) in S^. Projec- 
tion onto the dens ity, as in the interaction terms, then 
leads to Eqs. fl2.18| ), with 1/ti = Z7(q = 0), and l/r rei a 
weighted average over U (k— p). T re x is the single-particle 
relaxation or scattering time. 

Next we define a Grassmannian matrix field -B12, with 
1 = (m, a±), etc, by 



B 12 (x) =?7+(x)®?72(x) 



(2.19a) 



or, in Fourier space and with all components written out, 



? X™(q) = E^( k )^™( k + q) 



(2.19b) 



The subset of the Grassmann algebra that is given by the 
bilinear fields B is isomorphic to a set of classical, i.e. 
c-number valued, matrix fields. We exploit this isomor- 
phism, which in particular maps the adjoint operation in 
■0-space that we denoted by an overbar onto the complex 
conjugation operation, by introducing a classical matrix 
field Q, and constraining B to Q by means of a func- 
tional delta function.E3 We then use a functional integral 
representation of the latter, with an auxiliary or ghost 
field A that plays the role of a Lagrange multiplier, and 
integrate out the fermions. This way we obtain, 



D[ip, ip] e § ^ / D[Q] S[Q-B] 



= J D[ij\ e^W J D[Q] D[A] cxp (tr A(Q - B) 



D[Q]D[A] e 



A[Q,A] 



(2.20) 



In Eq. ( |2~20| ) we have defined an effective action, 



A[Q, A] = A dis [Q] + A mt [Q] + - Tr In ( Gq — iA 



ix tr (A(x) Q{x 



(2.21) 



Here and in what follows, tr is a trace over all discrete 
degrees of freedom that are not shown explicitly, while 
Tr is a trace over all degrees of freedom, including an 
integral over x. 



Go" 1 =-d T +V 2 /2m + f i 



(2.22) 



is the inverse of the free electron Green operator, and 
i t is clear from the structure of the Tr ln-term in Eq. 
(|2.21 ) that the physical interpretation of the ghost field 
is that of a self-energy. In writing Eq. ( 2.21 ), we have 
replaced the dummy integration variables A and Q by 
their transposeds, which carries a Jacobian of unity, and 
have anticipated the fact that Adis[Q T ] = Adis[Q], and 
the same for Aint- The latter two contributions to the 
action read 



Adis [Q] = Ajl2 [Q] + Af) [Q] 



-l 



2ttN F T! 



dx (trQ(x) 



1 



irNpTrel 



dx tr(Q(x)) ; 



and 



Amt [Q] — A int + A int + A 



» 



(2.23a) 



(2.23b) 



(2.23c) 



(2.24a) 



A 



E 

711,712,771 a 



[dxj:(-iy E 

J r=0,3 ni,n 2 , 
X [tr ((Tr®*o)OS"m+m(x)) 

x [tr ((r r ® S o)Q^ +m ,„ 2 (x)) 



A {t) - 



ttw 



dx e (-1)' 

r=0,3 



Til ,712 ,771 



x [tr ((T r ® Sl )Q^ ni+m {x)) 
x [tr ((T r ® Sl )Q^ +m ^(x)) 



dx E E 

r=l,2 ni ,ri2 ,m ex 



E 



x [tr ((r r « S „)Q™_ ni+m (x) 
X [tr ((r r ®S )<3-n 2) n 2 +m( x ) 



(2.24b) 



EE 



(2.24c) 



(2.24d) 



2. Properties of the Q-matrices 

We now derive some useful properties of the matrix 
field Q. Since B, Eqs. (2.19), is self-a djoint under the 
adjoint operation defined in Eq. (2.14b) and denoted by 
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a superscript '+', so is Q, Q + = Q. Notice that Q + 
is distinct from the hermitian conjugate of Q, which we 
write as = (Q T )* , with the star denoting complex 
conjugation. The definition of the former adjoint implies 
Q+ = C T Q T C, with 



ij /-tap _ r r 

^nm ~~ u nm o a fj t-i. 



and c defined in Eq. ( 2.14c ). We thus have 
Q = C T Q T C . 



(2.25) 



(2.26a) 



In addition, a direct calculation of the hermitian conju- 
gate of B, = B T , reveals another constraint that is 
inherited by the Q, namely 



-rgr- 1 



(2.26b) 



where the similarity transformation denoted by T has the 
property 



(rgr- 1 )„ m = o_ 



n — 1. — m — 1 



(2.26c) 



We now expand our mat rix fi elds in the spin-quaternion 
basis defined after Eqs. ( 2.14 ), 



Qi 2 (x) 
Aia(x) 



3 

E 

r,i— 
3 

E 

r.i— 



(r r ® Sj) rQi 2 ( x ) 



(r r ® *A 12 (x) 



(2.27a) 
(2.27b) 



where again 1 = ( ni, cti) , etc. In this basis, the relations 
expressed by Eq. ( 2.26a| ) imply the following symmetry 
properties, 



rQi2 — ( ) IQ21 j ( r — 0,3) , 

;Qi 2 = (-) r+l4 rQ 2 i . (r = 0,3; i = 1,2,3) 

, r <5 12 = r Q 2 i ' ( r = I7 2) , 

l r Q 12 = -lQ 21 , (r = l,2; 1 = 1,2,3) . 



(2.28a) 
(2.28b) 
(2.28c) 
(2.28d) 



Toge ther w ith the behavior under hermitian conjugation, 
Eq. ( 2.26h| ), this further implies 



r^l2 

l 0* 



"r 1 *— ni-1,— n 2 — 1 



aiQ 2 

' — Tlx — 1 j 



-712 — 1 



(r = 0,3) 
(r = l,2) 



(2.29a) 
(2.29b) 



Analogous relations hold for A by virtue of the linear 
coupling between Q and A. 

This completes the derivation of the action in terms of 
classical matrix fields. We reiterate that for disordered 
fermions, this action is adequate as long as one is inter- 
ested in physics dominated by soft modes, while for clean 
systems, a complete theory would have to keep an infi- 
nite set of interaction constants. We also recall that, by 
means of the Q-matrices, we keep explicitly only density 
modes in both the particle-hole and the particle-particle 



channels. While this is again physically well motivated 
in the disordered case, it is more problematic in the clean 
case. Due to the larger number of soft modes in the lat- 
ter, our procedure means that in certain clean system soft 
modes have been integrated out. While this integrating- 
out procedure is exact, it leads to undesirable properties 
of the effective field theory as we will discuss below. In 
the next subsection we will find further indications for it 
being advantageous to derive somewhat different effective 
theories for clean and disordered fermions, respectively. 



C. Saddle— point solutions and Gaussian 
approximation 

The physical degrees of freedom are now given by the 
matrix elements of the Q-matrices, and the physical cor- 
relation functions for number and spin density fluctua- 
tions can be expressed in terms of Q-corrclation func- 
tions. This implies that ultimately we will want to inte- 
grate out the ghost field A in Eq. ( [2.21 ). Before doing so, 
however, we examine a saddle-point solution to the field 
theory defined by Eqs. ( |2.2l[ ) - (2.24), and the Gaussian 



fluctuations about this saddle point. 



1. Fermi-liquid saddle-point solution 

It is possible to develop a systematic theory for the ex- 
pectation values of Q and A, rather than simply a saddle- 
point approximation for these quantities. We will come 
back to this point in Sec. Ill B 1| below. We further note 
that there are many saddle-point solutions, and that in 
general they have different symmetry properties. Phys- 



ically these different solutions correspond to states with 
different broken symmetries such as Fermi-liquid states, 
magnetically ordered states, or superconducting states. 
For now we concentrate on a-J^ermi-liquid state. In Sec. 
IIIB4 below, and elsewhere, c3 we will consider more ex- 



otic states. 

The saddle-point condition is 



SA 
SQ 



Qsp:A s 



SA 



(2.30) 



Q S p:A s 



The Fermi-liquid saddle-point solution is spatially uni- 
form, diagonal in frequency and replica space, and there 
are no nonzero expectation values that would describe a 
spontaneous magnetization, or a BCS-like order param- 
eter. In Fermi liquid-like phases, the saddle point values 
of both Q and A have the structures, 



j(x) = 612 S r o SiO Q n , 

sp 



;a 12 ( x ) 



S12 SrO <5i0 A„ , 



(2.31a) 
(2.31b) 
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and the average val ues o f Q and A have the same prop- 
erties. Using Eqs. ( 2.21 ) - ( 2.31 ), we obtain the saddle 
point equations 



i 

2V 



p 2 /2m + /i - iA n \ 



A,, 



nNpT, 



rcl 



(2.32b) 



Here e JWm0 is the usual convergence factor that re- 
solves the ambiguity presented by the equal-time Green 
function!] Notice that this saddle-point solution obeys 
all of the sy mmetr y properties of the Q that were de- 
rived in Sec. [I B 2 . For a physical interpretation of Eqs. 
( |2.32 ), we defin e a self-energy S„ = iA n . In terms of £, 
the Eqs. ( 2.32 ) can be rewritten, 



1 



1 



irNFT re i V 



E 

p 



\iu n - p 2 /2m + fj, - £„] 



l0 v E \- iuJm ~ p 2 / 2m 

p 

+/!-£,„r 1 . (2.33) 



This integral equation for the frequency dependent self- 
energy has a familia r str ucture. For vanishing interac- 
tion, r( s ) = 0, Eq. ( 2.33 ) reduces to the self-consistent 
Born approximation for the self-energy, or the inverse 
mean-free time, in a disordered fermion system. At zero 
disorder, on the other hand, Eq. ( |2.33| ) represents the 
Hart ree-F ock approximation for the self-energy. The full 
Eq. ( [2.33] ) we will refer to as the disordered Hartree- 
Fock approximation for the self-energy. It interpolates 
smoothly between the clean and dirty cases. 

In this context, a few remarks on ultraviolet divergen- 
cies and cutoffs are in ord er. In gener al di mensions, the 
sums or integrals in Eqs. ( 2.32 ) and ( 2.33 ) do not exist 
due to ultraviolet divergencies. Clearly, this is due to 
the long-wavelength and low-frequency approximations 
we have made by, for example, projecting onto density 
modes only. Any ultraviolet problems encountered in the 
theory are therefore artifacts, and should be cut off at 
some momentum scale k . In principle there are a num- 
ber of 'large' momentum scales that are candidates for 
ItQ, e.g. the Fermi wave number kp, the Thomas-Fermi 
wave number hrp, and the inverse elastic mean-free path 
l/l. For metallic densities, Uf ~ kpF, while kpl >> 1 
for weakly disordered systems, and kpl ~ 1 for strongly 
disordered ones. A popular choice is to use ko = 2tt/£ in 
any disopdered system, regardless of the strength of the 
disorder .t3 However, for weakly disordered systems this 
choice leads to qualitatively wrong answers for various 
nonuni-versal results, while using fco = kp is qualitatively 
correct SB For the universal phenomena that will be con- 
sidered in the present paper, on the other hand, the pre- 
cise value of ko is not important, and any of the above 
choices are acceptable. 



2. Gaussian approximation 

We now examine the Gaussian fluctuations about the 
saddle point discussed above. To this end, we write Q 
( p4[ ) as, 



(2.32a) an d A in Eqs. (|2~21 



A = A 



Sp 



SQ 
■SA 



(2.34a) 
(2.34b) 



and then expand to second or Gaussian order in the fluc- 
tuations SQ and SA. Denoting the constant saddle point 
contribution to the action by A sp , and the Gaussian ac- 
tion by Ag, we have 



A[Q,A] = A sp + A G [Q,A\ + AA 



(2.35) 



where A.A contains all terms of higher than quadratic 
order in the small quantities SQ and SA. The Gaussian 
action reads, 



Aq [Q, A] = A dls [SQ] + A mt [SQ] + ~ Tr (c sp SAG sp 



dx tr (<5A(x)5Q(x) 



SA 
(2.36a) 



with G sp the single particle Green function in saddle- 
point approximation, 



G sp (p, tu n ) = [iuj n - p 2 /2m + fj, - E n ] 



(2.36b) 



where £„ is given by the solution of Eq. (2.33). The 
spectrum of this saddle-point Green function determines 
the quantity Np, see Eq. (2.46) below. 



At this point it is important to notice that in the full 
actio n, Eq . ( p. 35 ) , the symmetry properties expressed by 
Eqs. ( 2.28| ) are enforced by means of delta- function con- 
straints. Truncating the theory at the Gaussian level 
spoils this property. Consequently, by using the con- 
straints in different ways before the truncation one ob- 
tains different Gaussian theories. Here we opt to use the 
constraints to rewrite the theory entirely in terms of ma- 
trix elements SQ12 with n\ > n%. Correlation functions of 
matrix elements that do not obey these frequency restric- 
tions are related to Eqs. (2.37) below by means of Eqs. 
(Uf). As we w 11 see, this choice results in a Gaussian 
theory that reproduces results that are well-known from 
man y-body perturbation theory. With this procedure, 
Eq. ( 2.36a ) is a quadratic form that can be diagonalized 
using the techniques discussed in Ref. |l6|. One obtains 

\ V 

•W)l 2 (Pl)iW) 3 4(P2)/ G = ,-P2 $rs Sij 

x;M^ 34 ( Pl ) , (2.37a) 



;(M) 12 ( Pl )i(<5A) 34 ( P2 ))^ = - ( 5 Pl 



V 



-P2 ^rs $ij 



X r N J3.M ■ ( 2 ' 37b ) 



9 



denotes an average with a weight Ag- In _ p 



Here L_ 

Eq. ( 2.371: ), we have defined a correlation function for 
the fluctuations of a field, A, that is decoupled from the 
Q-correlation function. A is defined by 



A12 = - f ni n 2 Al2 + Ql2 



(2.38a) 



with 



wm (k) = — G sp (p, bj n ) G sp (p + k,cj m ) . (2.38b) 



The matrix inverse M 1 for the particle-hole channel is 
given by, 



3-4 0,3 



<5l3 ^nmj (p) 



^Q 1 a 2 <5q3Q4 ' X? ninj (p) 2?n 3 n4 (p) 

4(5 r0 (5. t o 



0) 



(2.39a) 



with 



* i3 / 12 = 1 + <5 til „ 2 (1 - 2S r3 ) (Jio - 51 . ( 2 - 39b ) 

and r<°) = -r( s ), and rt 1 ' 2 ' 3 ) = rM. For the particle- 
particle channel, one finds, 



1,2-^12*34 ( P) = - (5 l+2,3+4 5l3 1,2^2 ^lU (p) 

- ««j <W 2 ^ 3 «4 4Tr( c )pW n2 (p) P„ 3 „ 4 (p) 



where 



1,2-^2 - 1 + ^nin 2 (<^o - 

j=i 



(2.39c) 



(2.39d) 



Similarly, one finds for the matrix inverse N 1 in the 
particle-hole channel, 



0,3-^12,34^) = _<5 13^24 0,3^2 ¥>nm 2 (p) 

and in the particle-particle channel, 

1,2-^12,34 (p) = ^13 <*24 1,2 J 12 ^«in 2 (P) 

In the preceding equations, 

1 



(2.39e) 



(2.39f) 



and, 



1 



TTNpTr 



Vmm (p) 



(2.40a) 



1 + 2TT (i) 

^ ] f n,n+n 2 — "i (P; 



(2.40b) 



with X>(°) = 2?( s \ and £>W = X>( 2 > = X>( 3 ) = X>W and 



P&Up) (p) 



1 + 2Tr (c) 

^ ' 2^n,ni+ri3— n(p) 



(2.40c) 



D. Physical correlation functions 

We conclude this section by interpreting the results ob- 
tained in the preceding subsection, and expressing some 
physical quantities of intere st in te rms of Q-correlations. 

First of all, ip nm (]&), Eq. ( 2.38b ), is related to a disor- 
dered Lindhard function. We write 



i^nm(k) = Q(nm)$ + ®(—nm)ip(k, 0„- m ) 



(2.41) 



with Q n — 2itTn a bosonic Matsubara frequency. For 
both clean and disordered systems, <E> approaches a fi- 
nite constant in the limit of vanishing wave number and 
frequency. On the other hand, the behavior of ip in the 
same limit depends qualitatively on the presence or ab- 
sence of disorder. For dirty systems, to leading order as 

T re l — ► OO, 



ip{\k\ -> 0,fl n -» 0) = irN F T re l 



(2.42) 



while for clean systems, ip is singular in the long- 
wavelength, low-frequency limit. For example, in d = 3, 

/Hi no r>\ i7rNF rn U i^n + \k\v F 
tp(\k\ -> 0,O„ -> 0) = — - — sgn(S2„)ln T 



2|k|t>F 



fir, 



\k\v F 
(2.43) 



Comparing Eqs. (2.42) and (2.43), we see that the l/|k 



and 1/Cl n singularities that occur in the clean case are 
cut off by a mean-free path, and a relaxation time, re- 
spectively, in the disordered case. This implies that the 
A- A correlation function is soft for clean systems, and 
massive for dirty ones. This will turn out to be a funda- 
mental difference between clean and dirty systems, which 
reflects the fact that single particle excitations are soft 
in a clean system. These observations indicate that for 
the further development of the theory it will be useful to 
treat clean and disordered systems in somewhat different 

ways. 

ip determines the propagator T>, Eq. ( |2.40a ). For finite 
disorder, n\ri2 < 0, and in the long-wavelength, low- 
frequency limit one finds 
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f ni n 2 (p) 



Dp 2 + |n ni . 



(2.44) 



with D VpT re i/d a bare or Boltzman n diff usion 



coeffici ent .Ej The propagators given by Eqs. (2.37a) and 
( |2.39aD for n x n 2 < are then the soft or massless prop- 
agators of the nonlinear sigma-model that is used to de- 
scribe the metal-insulatoiU,ia|nsition in disordered inter- 
acting electronic systems Ota 

Second, the physical correlation functions of the 
single-particle spectral density, the particle number den- 
sity, and the spin density, can be easily expressed in terms 
of the Q-correlation functions. By keeping an appropri- 
ate source in the action while performing the transfor- 
mation to Q-variables, we obtain an expression for the 
density of states, N, as a function of energy or frequency 



N(e F +uj) = ~-Re(° Q nn (x 



lUJn — >LJ-\~iQ 



(2.45) 



In saddle point approximation, we have for the density 
of states at the Fermi level, 



N F = 



-2 1 
~ V 



E Im G sp (p, iuj n -> iO) 



(2.46) 



This defines the density of states that has been used 
throughout as an energy normalization factor. Similarly, 
one finds that the number density susceptibility, \ n , and 
the spin density susceptibility, Xsi are given by, 



X «(q,r! n ) = 16T 



E E (t(AQC_„, mi (q) 

mi ,?7i2 r— 0,3 

x ;(AQC 2 , m2+ „(-q)) , (2.47) 



with AQ = Q - (Q) , and x 



(0) 



Xn and x 



(1,2,3) 



Xs 



By substituting the Gaussian propagator we see that the 
Gaussian theory yields x„ and Xs m & disordered RPA- 
like approximation. Notice that, for w eak dis order, the 
last term on the right-hand side of Eq. ( 2.39a ) is smaller 
than the first two terms by at least a factor of 1/ri, 
and therefore does not appear in the usual resummation 
schemes of many-body perturbation theory. Neglecting 
it, one finds, for small |q| and f2„, 



x«( q ,fi n ) = x« 



|n n |+DWq 2 



where 



xi; } =iv F /(i-7v F r 



(2.48a) 



(2.48b) 



is the static susceptibility in the long-wavelength limit, 
renormalized by the interaction, and 



L> w = D (1 - N F T {l) ) 



(2.48c) 



is the renormalized diffusion coefficient. Here we have 
used the identity T^ n D ran (q — * 0) = —Np/2, which 



expresses partic le nu mber conservation. Notice that for 
i = 1,2,3, Eqs. ( 2.48 ) describe a disordered Stoner crite- 



rion: For NfTM = 1, the spin diffusion coefficient van- 
ishes, and the static spin susceptibility diverges. In the 
clean limit one recovers RPA proper, 



X (i) (q,fi„) = 2 X o(q,fi„) [l + r«2 Xo (q, fi„) 



(2.49) 



where Xo(q, H n ) = TJ2 m ¥'m-n,m(q) is the free electron 
density susceptibility per spin. 



III. THE DISORDERED FERMION SYSTEM 

In this section we concentrate on the disordered sys- 
tem, and in particular perform a symmetry analysis in 
order to identify the soft modes in our system. This will 
allow us to explicitly separate the soft modes from the 
massive ones, and to formulate an effective field theory 
for which the soft modes remain manifestly soft to all or- 
ders in perturbation theory. In Appendix |bJ we point out 
that all of the considerations and results of this section 
have precise analogies within the ^-representation of the 
O(N) Heisenberg model. These analogies are often help- 
ful, since ^> 4 -theory is technically much simpler than the 
matrix field theory under consideration here, and hence 
it displays the basic structure of the symmetry analysis 
and its consequences more clearly. 



A. Symmetry analysis 

1. Basic transformation properties 

Let us perform a symmetry analysis of our field theory. 
To this end, we start with the action in fermionic form, 
formulated in terms of bispinors, Eqs. ( 2.16Q - (2. IS). Our 



plan of attack is to consider the various terms in the ac- 
tion separately. The crucial symmetries will involve the 
free electron part, So, and the d isorder part, Sdis- The 
interaction part, Si n t, Eqs. ( 2.17 ), will turn out to always 
be effectively proportional to a frequency, at least near 
the disordered Fermi-liquid FP that we will identify be- 
low. Since we are interested in low-frequency effects, this 
will ultimately imply that Si n t does not change the con- 
clusions with respect to which two-particle excitations 
are soft or massless that are reached by considering the 

noninteracting action. 

For zero external frequency, iuj n — in Eq. (2.16a), 
the free fermion action is invariant under transformations 
that leave invariant the expression 



EE tr (V)s®<) = M) 



(3.1) 



ii 



with (77 , 77) a generalization of the scalar product defined 
in Eq. ( 2.15| ) that includes summation over the frequency 
and replica indices. Fro m the structure of the disorder 
part of the action, Eq. ( 2.16b ), it is clear that Sdis is 
also invariant under transformations that leave invariant 

(v,v)- 

Now let 77 be transformed by means of an operator T: 
i] — > Tr). Remembering that r] + is relate d to 77 by means 
of the charge conjugation matrix G, Eq. fl2.25| ), and using 
C T = G -1 , we find that in order to leave (77, 77) invariant, 
T must obey 



T CT = C 



(3.2) 



For a system with 2N frequency labels (N positive ones, 
i nclu ding 0, and N negative ones), and n replicas, Eq. 
( |3.2|) defines a representation of the syaiplectic group 
Sp(8Nn,C) over the complex numbers C.Eil Throughout 
this section we will denote the real, complex, and quater- 
nion number fields by 1Z, C, and Q, respectively. In what 
follows we will need to consider only small subgroups of 
the 8Nn(8Nn+ l)-parameter group Sp(8Nn,C). 

Now we ask what transformations of our composite 
variables Q correspond to the above transformations of 
the bispinors. Since the Q-matrix is cons train ed to the 
bilinear fermion field B = rj + (8)77, see Eq. ( 2.2C ), 77 — > TV? 
implies Q -> TQT-\ with T = CTC T . The ghost field 
A couples linearly to Q and B, and therefore transforms 
like Q, i.e. by means of the matrices T. T oget her with the 
T, the transformations T also obey Eq. (3.2), and there- 
fore also form a representation of the group Sp(8iVn, C). 

Now we consider a specific transformation, namely a 
rotation in frequency space given by, 



T nm = S i0 S r0 S a /3 \ Km 1 + (<W + C 2 ) (COS 0-1) 



+ (5nn 1 S n 



SiQ S r Q 5 a @ i Tl 



(3.3) 



This transforms a pair of spinors with frequency labels 
ni and 712, respectively, into linear combinations of the 
same pair, wit h a mixing angle 9. These transformations 
obey Eq. ( |3.2| ), and thus are elements of Sp(8Nn,C). For 
fixed ni and n-i they represent an SO (2) subgroup of 
Sp(8A r 77, C). The corresponding transformations T of the 
Q and A-matrices are identical with the T, T = T. Under 
an infinitesimal clement of this subgroup, the Q-matrices 
transform like 



(3.4a) 



with 



Qnm)0 + O{6 2 ) . (3.4b) 



Here we have shown only the frequency indices, since all 
other degrees of freedom are unaffected by the transfor- 
mation. The A-matrices transform accordingly. 

O f the various pieces of the action in Q-A formulation, 
Eq. (2.21), Adis and / dx. tr (AQ) are invariant under the 
above transformation, but Tr h^G^" 1 — iA) and Ai n t are 
not. For now, let us consider the transformation of the 
former: 



Tr In G, 



iA 



Tr In G 



iA 



'Tr(Gfe) , 
(3.5a) 



with G = (Gq 1 - iA)- 



and 



(3.5b) 

The action of the noninteracting Fermi system therefore 
transforms like 



Ao — A Ai' 



An + SA = A + - Tr (G Siu) 



(3.6) 



We will con sider the transformation properties of Ai n t in 
Sec. IIIIA3I below. 



2. A Ward identity for noninteracting electrons 

Because of the complexity of the interacting case, and 
for the reasons noted at the beginning of the previous 
subsection, we first discuss the noninteracting system. 
Let us introduce a source J for the Q-fields, and consider 
the partition function 



Z[J] 



D 



[Q] D[A] e Ao+ I dx tr ( J(x) Q(x) ) . (3.7) 



By performing an infinitesimal transformation, Eq. 

on the Q and A-fields, one obtains from Eq. ( |3.7| ) 
the identity, 



0= / D[Q]D[A] 



SA + j dx tr (J(x) 5Q(x)) 

.0+ J dx tr (,7(x) Q(x)) ^ 



xe 



Differentiating this identity with respect to (jJ^" ni and 
putting J = yields 



= (5Ao° Qf in2 (x 



\(SQ)ZnM 



(3.8b) 



where (. . .) denotes an average by means of the action 
Ao- From the explicit expression for SAq, Eq. (|3~6|), we 
see that the first average is essentially a frequency times 
(GQ). By using the identity^ (GQ) = -2i(Q 2 ) we fi- 
nally obtain 
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(3.9) 



This is. t>ha -desired Ward identity for noninteracting 
electrons ,EjO'E£| which relates a two-point Q-correlation 
function to the one-point function. For sgnni = sgnn2, 
the right hand side vanishes as rt\ — n 2 —> 0, and 
the identity tells us nothing inter estin g. However, for 
sgnni ^ sgnr«2 we see from Eq. (2.45) that the right- 
hand side approaches TrN(ep)/2, which is nonzero every- 
where inside the band. The correlation function on the 
left-hand side must therefore diverge. Due to rotational 
invariance in replica space we have (Q a/3 ) ~ S a p, and 
(Q al3 Q~* 5 ) ~ 5 ai 5f} S . We finally obtain, for nm 2 < 



and |O ni . 



0. 



k=0 



7rN(e F ) 

I6\fl ni -n 2 



(3.10) 



Here N(cf) is the exact density of states at the Fermi 
level of the noninteracting electron system. The salient 
point is that, as long as N(cf) > 0, the Q-Q correla- 
tion function at zero momentum diverges like l/|f2 ni _„ 2 1. 
Equations ( f2.37a[ ), fl2.39a| ), and ( |2.44j ) show that the 
Gaussian propagator has this property. The Ward iden- 
tity ensures that it holds to all orders in perturbation 
theory. We have therefore identified oQ""„ 2 for n\ri2 < 
as a soft mode. Rotational symmetry in replica space im- 
plies that the oQ"m' nm < 0; a A are & l so s °ft- This 
can also be seen by applying the above rotation procedure 
in frequency space directly to spinor pairs with different 

replica indices. 

From the Gaussian propagator, Eqs. ( [2.37a ), (2.39), 
as well as from physical considerations, one expects all 
channels to be soft, not just the r = 0,i = channel. 
This is indeed the case. These additional soft modes are 
not controlled by sepa rate Ward identities, but rather 
are related to Eq. ( |3.9| ) by additional symmetries of the 
action. In order to see this, we consider transformations 



T, 



[Snn 2 X r (T r <g> Si) + (1 - S nn2 ) (t <g> S )] , 

(3.11a) 



with xo = Xa = 1, and x\ = x 2 = i, and n 2 some fixed 



frequency index. These T obey Eq. (3.2), and in addition 
they are unitary. It is easy to check that under these 
transformations Q transforms like 



(n ^ n 2 ) 



(3.11b) 



If the action is invariant under such a transformation 
for given i,r, then (IQIQ) = (qQoQ), and hence the 
two-point Q-correlations in the i, r channel are also soft. 
Obviously, the noninteracting actions is indeed invariant 
under these transformations, so the two-point correla- 
tions arc soft in all channels. 



Note that Eqs. (3.11) imply that the spin-singlet (i = 
0) and spin-triplet (i — 1,2,3) channels are equal, and 
so are the particle-hole (r = 0,3) and particle-particle 
(r = 1,2) channels. In particular, they provide a techni- 
cal explanation for why the so-called Cooperons, i.e. the 
two-point Q-correlation functions in the particle-particle 
channel, are soft modes. We also mention that in the 
clean case, there still is only one basic Ward identity, 
but many moi£ additional symmetries than in a disor- 
dered systemO These additional symmetries result in 
many additional soft modes. This feature leads to the 
somewhat paradoxical fact that it is easier to derive an 
effective field theory for disordered electron systems than 
for clean ones. 



3. A Ward identity for interacting electrons 

Now we add the electron-electron interaction to our 
considerations. For simplicity, we restrict ourselves to 
a discussion of the particle-hole spin-singlet interaction; 
the discussion of the remaining interaction channels pro- 
ceeds analogously. 

Let us consider again an infinit esim al element of the ro - 
tatio ns in frequency space, Eq. ( |3.3| ). From Eqs. (2.241) 
and ( 3.4b ) we find that the particle-hole spin-singlet part 



of the action transforms like A 



A%1 + 5 A£>, with 



» 



7 r=0,3 n'n' 



O T7 M - u O" M 

r^m + ttij-n;)^ \ ' r^ni,ri 2 -(n' 1 -n 2 ) W 



i rQn2 + (n' 1 -n' 2 )ji 1 ( X ) + ^n 2 ,ni-(n;-n;)( X ) 



(3.12) 



Now we follow the same steps that led to Eq. (3.81), 



except that we differentiate with respect to o jf"„ 3 with 
n3 > 0, ri4 < 0. This yields 



+ («5A„t8Qi 4 (x)) . (3.13) 



Choosi ng ni > 0, n 2 < 0, and using (Q nm ) ~ 5 nm as well 
as Eq. ( |3.12 ), we obtain the Ward identity in the form, 



W mt + 80 ni _ n2 £ fdy (°Q^ n2 (y) Stef„ 4 (x) 

7 J 

= S nin3 5 n2ni ^oQnini( X )) — (oQn 2 n 2 ( X ) )) 



(3.14a) 



where, 

W int 
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1"2 



O 77 (V) 



— °n 77 ( \ _1_ 0/)T7 / n 

— r0n 2 ,n 1 -(ri' 1 -Ti2)(y) 

(3.14b) 

with m, ri3 > 0, and 7x2, Ti4 < 0. 

We now ne ed to analyze the 3-point correlation func- 
tions in Eq. (3.14b), and ascertain that the Q ni n 2 with 
niri2 < remain soft modes in their presence. To this 
end we first split up the Q into their averages, and fluctu- 
ations: Qnm — (Qnm) ~~ (^Q)n-m — ^nm(Qnn) ~t~ (^Q)nm- 

Since ni 7^ 77,2, and 77,3 7^ 714, we have (Q ni n 2 ) — 
(Qn 3 ni) = 0. Furthermore, if we put nl-y = n' 2 in the 
3-point fun ctions, then the expression in square brack- 
ets in Eq. ( 3.14b ) vanis hes, s o effectively we also have 
(Q n ' in ' 2 ) = 0. Equation (3.14a) then takes the form, 



W mt + 8Q ni _ n2 £ J dy (%(AQ)l] n2 (y) 



, a/3 



xg(AQ)^(x)^ 

= <W 3 <W 4 ((o< ni (*)) - (g0i 2 (X))) , 

(3.14c) 

and Wi n t can be written, 

w mt = -lerw^ ((°q- i(x) ) - (°q- 2 ( X ))) 

7 

^ ' (^n'j— n' 3 ,n,2— ni 3n' 1 —n 2 ,n 1 — n 2 ) 
n 'l n ' 2 

x /dy(g(AQ) 7 7 n ,(y)°(AQ)i 4 (x)) 
-16TWW«rE T E 

7 r=0,3 njiij 

x/g(AQ)ifx)°(AQ)I7 n ,(y) 



°(aq)i: 



n 2 - (n'j— Ti^ ) 
77 



(y) 



-?(AQ)£ +K _ n , )ina (y) 



-°(AQ)I 
-?(AQ); 



2 +(n^— n 2 ) 



,m(y) 



2,ni-(n' 1 -ri2) 



^(y) 



(3.14d) 



Next we remember that the composite variables Q cor- 
respond to products of fermions variables, e.g. Q ~ ipip 
in the particle-hole channel. Specifically, 



oQriin 2 ~ „ ^ ] {' l Pni,(7' l Pn 2 ,(T + VV12 ,<r VVii ,ct) 



(3.15a) 



(3.15b) 



Therefore, writing the correlation functions in terms of 
AQ ~ V'V ~ (V'V') makes them Q'hreducible, but not 
^-irred ucibl e . For instance, the 2-point function in Eqs. 
( |3.14c ) and (3.14d) has the structure, 



E 7 / dy (°(AQ)I 7 „ 2 (y) g(AQ)^„ 4 (x) 

— A a f\ f\ X^ -A- t) X^ 

(3.16) 

The replica structure of this identity is such that X^ 1 ' 
and X^ are replica independent £a We can therefore 
suppress the replica index in what follows, as we have 
done al ready in Eqs. (3.15). The frequency st ructur e 



of Eq. ( 3.16 ) is very general: By virtue of Eq. ( 3.15a ), 
the full 2-point function is automatically proportional 
to 

Sni-n 2 ,n 3 -n4, due to time translational invariance, 
and the disconnected part is in a ddition proportional 
to S ni n 3 S n2 n 4 ■ Using Eqs. ( 3.15 ) in the 3-point Q- 
correlations shows that they also contain parts that are 
proportional to S nin3 S n2ni , and others that are propor- 
tional to <5ru-n3.n2-TJ4- This suggests to break up the 
Ward identity into two separate pieces with these two 
frequency structures, respectively. 

By collecting the above arguments, we find that the 
coefficients of 6 nin3 S n2n4 in the Ward identity obey the 
relation, 



oQmm ( x ) 



^QtnM)) = ^[^ ni ^ n2 X^ n2 

-N F {N F T^fY^J . (3.17a) 



Here Y W is related to the piece of the 3-point correlation 
functions that is proportional to S nin3 S n2n4 . An explicit 
calculation yields, 



y(i) = ~ 8 

1 n 1 n 2 N 3 T {s) 



x <j G(x, y; w ni ) ( (ip^ (y)s ip n ' 2 (y)) 

X ((■0T l2 -n' 1 +n 2 (y)so0r l2 (x)) 



-G(x, y; uj n2 ) (?/>„/ (y)s Vv» 2 (y)) 

X (('0n 1 -» l ' 1 +n 2 (y)soV , n 1 (x)) 



dis 

(3.17b) 

Here the (ipipip?p} c are '(/'-connected correlation functions, 
with the cumulant taken with respect to the quantum 
mechanical expectation value only, and G(x,y;cj„) = 
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(■0n,o-(x) ipn,a(y)) is the Green function for a given disor- 
der configuration. For Tm — > +0 and Tn-i —* — 0, so that 



by means of additional symmetries that are not broken. 
First of all, rotational symmetry in replica space again 
0, the left-hand side of Eq. ( |3.17a| ) approaches implies that the qQ"^, nm < 0, are also soft.0 Next we 

show that the spin-triplet channel also remains soft. To 
this end, let us denote the (identical) coupling constants 
in the three branches of the spin-triplet channel by f\ , 
i = 1,2 , 3. Now we pick from the transformations T, Eq. 
d3.11a| ), the one that interchanges, say, i = and i = 1. 
The action is still invariant under this transformation, 



^Tii — n\ 

again 7rAf(e F )/2, with -/V(e^) the single-particle density 
of states at the Fermi level. However, since we are now 
dealing with interacting electrons, N{tp) includes inter- 
action as well as disorder effects. Ynl n2 in that limit 
approaches some number ttY^ /2. We then obtain, for 
small fL,_„,, 



x<n = — 

ni " 2 i6|n ni . 



N(e F ) + n f (n f t^) 2 y^ 



(3.18) 



Notice that the -0-connected correlation functions vanish 
for noninteracting electrons, so that in an expansion in 
powers of the interaction constant, is of 0(1). The 
term in brackets on the right-hand side of Eq. ( 3. 18] ) 

is equal to 7r(7V(e F )) 2 /8iJ , wiltLjEf the quasi-particle or 
specific heat density of states.B'HL3 

Similarly, one finds for the coefficients of <5 ni _„ 3in2 _„ 4 
in the Ward identity, 

"ni-12 A nin 2 ,n 3 n 4 — iVj?^iVi?l ) * ni n 2 , 113114, 

-16T^T { [G(x, x; w Wl ) - G(x, x; to n2 )] 

(x,xH , (3.19a) 

J dis 

where 

(j) n3n4 (x,x) = J dyG{y,x;uj n3 )G(x,y;uj ni ) , (3.19b) 

and Y ^ is again given in terms of connected 4rip correla- 
tion functions. If desired, it can be calculate d in p ertur- 
bation theory in r( s ). For tt niri2 -> 0, Eqs. (|3.19|) yield 



y (2) 



T/\n ri 



Again , the G a ussia n 



propa gator has that property, see Eqs. ( [2.37a ), ( 2.39 ), 

&■ . 

From Eqs. ( 3.18 ) and ( 3.1 9| ) we see that gQ"^ remains 

a soft mode in the presence of interactions. This is due 
to the frequency structure of the interaction term, which 
ensures that the 2-point Q-correlation function remains 
soft rather than acquiring a mass, as one might naively 
expec t from an inspection of the action, Eqs. (2.21) - 
( [2.24 ). Note that this softness is not restricted to the 
particular linear combination of Q that constitutes the 
particle number density. It therefore is not related to 
the particle number conservation law. Rather, it is a 
consequence of the spontaneous breaking of a continu- 
ous symmetry in the system, viz. the invariance, at zero 
external frequency, under the rotations between positive 
and negative frequencies discussed above, or between re- 
tarded and advanced Green functions. The soft modes 
in question are the corresponding Goldstone modes. 

As in the noninteracting case, there is a variety of other 
soft modes that are related to the above Ward identity 



provided that we also interchange and 1^ . This 



shows t hat \ Q a lso ob eys Eq. (3.16), but with the 
in Eqs. fl3.17| ) - fl3.19a| ) replaced by rW. Since this can 
be done for all branches of the spin-triplet channel, it 
follows that they all are soft modes. Alternatively, we 
can invoke rotational invariance in spin space (an SU (2) 
subgroup of the large symplectic group), which implies 
that all three components of the spin-triplet have identi- 
cal correlation functions. This argument also shows that 
there are no singlet-triplet cross correlations. By means 
of similar arguments one easily convinces oneself that the 
particle-particle channel is still soft, with the appropriate 
Cooper channel interaction ampli tude re placin g or 
rW. For the action shown in Eq. fl2.2l|) - fl2.24| ), we thus 
conclude that all channels are still soft, and all t he sof t 
modes can be traced to the Ward identity, Eq. ( 3.14a ), 
and additional symmetry relations. If some of the addi- 
tional symmetries are broken externally, e.g. by means 
of an external magnetic field, or by magnetic impurities, 
then the respective modes acquire a mass, see Refs. [24 16 



4- Separation of soft and massive modes 

From the previous subsection we know that the corre- 
lation functions of the Q nm with nm < are soft, while 
those with nm > are massive. Our next goal is to 
separate these degrees if freedom in such a way that the 
soft modes remain manifestly soft to all orders in per- 
turbation theory. To give a simple example, in an O(N) 
4> -theory this is achieved by writing the A^-component 
vector field as (f>(x) — p(x) 4>(x), with p the norm of the 
vector, and |</>(x)| = 1. With this separation, only gra- 
dients of <f> a ppe ar in the theory, and so <p is manifestly 
soft, see Ref. [35] and Appendix |b]. We are looking for an 
analogous separation of our matrix degrees of freedom. 
For noninteracting electrons this has been done in Refs. 
|l7| and [l^, and we will construct a suitable generalization 
of the procedure used by these authors. 

The matrices Q that we are con sideri ng are complex 
8Nnx8Nn matrices that obey Eqs. (2.26). Alternatively, 
we can consider the Q as 2Nn x 2Nn matrices whose 
elements are spin-quaternions, or 4Nn x 4Nn matrices 
whose elements are quaternions. For our purposes the 
last choice will be most convenient. We thus study a set 
of matrices Q with quaternion valued matrix elements 



' nm,c7(7' 



(n,m — 1,...,N; a, (3 = 1, 



,n; a, a =| 
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Q 



a0 



r=0 



that obey 



Q = C 1 Q 1 C 



(3.20a) 



(3.20b) 
(3.20c) 



with r from Eq. ( [2.26c ), and C = 1<8>T2, where 1 denotes 
the ANn x 4^77, unit matrix. This set of matrices is 
invariant under transformations T that-ebey TCT T = C, 
i.e. under symplectic transformations £3 

For what follows, the constraint expressed by Eq. 
(3.20c) is somewhat awkward to handle. We therefore 
analytically continue the matrix elements of Q, consid- 
ered as functions of the imaginary frequencies, iu) n and 
ito m , to the complex plane: iuj n — > z\, and iui m — > z%. As 
functions of Z\ or z%, the Q have a branch cut on the real 



axis, and from Sec. [II B 3 we know that those matrix 



elements with Z\ and zi approaching the real axis from 
opposite sides are soft modes. The analytic continuation 
of Eq. Q3 20cD reads 



(Q^)ziZ 2 — Qz*Z* 



(3.21a) 



In particular, for real frequencies, u> and a/, we have 



(Q^)w±iO,w'±iO — —Qui^iQ.u'^iQ 



(3.21b) 



If we analytically continue Q onto the unphysical Ric- 
mann sheet, Q UtU > is thus formally anti-hermitian. It 
therefore has ANn imaginary eigenvalues Aj-to (Xj <E 
ilZ; j — 1, . ..,4jVn), and can be diagonalized by 
means of unitary transformations S E U(ANn, Q) (i.e. 
by unitary ANn x ANn matrices whose elements are 
quaternions). Now U(ANn, Q) is isomorphic to the 
unitary symplectic group USp(8iVn,C) = \J(8Nn,C) n 
Sp(8A r n, C) .EiJ This means that the Q can be diagonal- 
ized by means of unitary matrices that are also symplec- 
tic, and hence leave the set of Q invariant. That is, the 
most general Q can be written 



Q = S D S' 



(3.22) 



where D is diagonal, and S £ USp(8A^n, C). 

However, diagonalization is more than we want. Since 
we know that the Q nm with nm < are soft, while those 
with nm > are massive, we are interested in generat- 
ing the most general Q from a matrix P that is block- 
diagonal in Matsubara frequency space, 



P : 



P> 
P< 



(3.23) 



where P > and P < are matrices with elements P nm where 
n, m > and n, m < 0, respectively. This can easily be 
achieved. Since the analytic continuations of P y and P< 



are anti-hermitian, the most general P can be obtained 
from D by an element U of USp(4iVn, C) x USp(47Vn, C). 
The most general Q can therefore be written 



Q = S P S~ 



(3.24) 



with S = SU . The set of transformations S is 
the set of all cosets of USp(8JVn, C) with respect to 
USp(47Vn, C) x USp(4iVn, C), i.e. the S form the homoge- 
neous space USp(87Vn,C)/USp(4iVn,C) x USp(47Vn,C). 
The corresponding most general Q on the imaginary fre- 
quency axis is generated from the most general P by a set 
of transformations that is isomorphic to this coset space, 
and that can be explicitly constructed by reversing the 
Wick rotation that led us from imaginary frequencies to 
real ones. For our purposes, we will not need this ex- 
plicit construction, and we will not distinguish between 
the two isomorphic spaces. 

This achieves the desired separation of our degrees 
of freedom into soft and massive ones. The massive 
degrees of freedom are represented by the matrix P, 
while the soft ones are represented by the transforma- 
tions S e USp(87Vn,C)/USp(47Vn,C) x \JSp(ANn,C). 
To come back to the O(N) example at the beginning 
of this subsection, the analogy is as follows. For the 
iV-component vector, the direction is a point on the 
(N — l)-sphere, which is isomorphic to the homogeneous 
space 0(N)/0(N - 1) X 0(1), see Appendix g. The 
unitary-symplectic coset space identified above is a ma- 
trix generalization of this. 

In order to formulate the field theory in terms of the 
soft and massive modes, one also needs the invariant mea- 
sure I[P], or the Jacobian of the transformation from the 
Q to the P and the S, defined by 



Jd[Q]... = Jd[P]I[P] Jd\s] 



(3.25) 



Since we will not need to know the explicit form of I[P] 
in what follows, we relegate its derivation to Appendix 



5. A Ward identity based on a local symmetry 

In addition to the Ward ide ntity de rive d on th e basis 
of a global symmetry in Sees. Ill A 2 and III A 3 above, 
for the derivation of an effective field theory in the next 
subsection we will also need an identity that is based on 
a local symmetry. Such relations with the structure of a 
Ward identity or a Noether theorem can be derived on 
the basis of either the replica or the gauge structure of 
the theory. For our present purposes, we consider a lo- 
cal U(l) gauge transformation of the Grassmannian field 
theory, Eqs. (2.2), that is independent of imaginary time 



-£a(x) 



(3.26) 
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We define a (d+l)-component current vector q u = (go, q) 
with 'spatial' components (i = 1, . . . d) as 



9o O) = ^^(aOVvOz) , 



(3.27a) 



q(x) 



-1 
2m 



(i(i)Vi(i)+i(i)Vi(i)) 



(3.27b) 



and add to the action a term Sa that describes a coupling 
of the current to a r-independent vector potential A u — 
(A ,A), 

S A = ~i J dxA^)q^{x) - / dxA 2 (x)g (x) . 

(3.28) 

It is then straightforward to establish that the partition 
function is gauge invariant, i.e. it is invar iant u nder the 
local gauge tranformation given by Eq. ( |3~26| ), supple- 
mented by a transformation of the vector potential, 

A„(x) ^ A„(x) + d M a(x) , (3.29) 

with <9 M = (d T , V). The invariance statement is 



Z[A U ] = Z[A u + d u a] 



(3.30) 



A Taylor expansion of Eq. ( 3.3C ) immediately leads to 
the conservation law 



(3.31) 



d tl [{q tl (x))--A^x){q (x))\ = 



where (...) denotes an average with respect to the ac- 
tion S + Sa- By putting A n = 0, one simply obtains the 
continuity equation for the fermion number density. An- 
other useful identity can be obtained by integrating Eq. 
(3.31) over space and time, differentiating with respect 
to Ai(y), and then setting A n = 0. The result is, 



d^„ j(7 (x) 



V>n,"-( x ) 



^m,a' (y) 



(Vv<x( x ) ^„ )<T (x)) 



(3.32) 



This identity holds for a particular realization of the dis- 
order. If we perform the ensemble average it still holds, 
with the brackets now including the disorder average in 
addition to the quantum mechanical expectation value. 
One then has a relation between a homogeneous current- 
current correlation function and the particle number den- 
sity that is known as the f-sum rule. The above deriva- 
tion makes it clear that it is closely related to the particle 
number conservation law. Within models or approxima- 
tions that describe an effective single-particle problem, 



the qua ntum mechanical average on the left -han d side 
of Eq. (H|) factorizes. Furthermore, Eq. ( |3~3^ ) then 
holds for each term in the n-summation separately, since 
in a single-particle problem there is no frequency mix- 
ing. Within approximations where the disorder average 
factorizes as well, the identity can be written in terms of 
disorder averaged one-particle Green functions as, 



SijG(x 



0,w„ 



— [ dy [diG(x - y, w n )] 
to J 

x [djG(y - x,u/„)] 



or, in Fourier space, 



q 



(3.33a) 



(3.33b) 



This identity will be useful later. 



B. Effective field theory for disordered fermions 

We are now in a position to formulate the theory in 
such a way that the soft and massive modes, respectively, 
are separated, and that the former remain manifestly 
soft to all orders in perturbation theory. This formu- 
lation also provides a technically satisfactory derivation 
of the generalized nonlinear sigma-model for interacting, 
disordered fermions. We will first discuss this deriva- 
tion, as well as corrections to the sigma-model. Then we 
will discuss the Fermi-liquid FP and the critical FP that 
describes the Anderson-Mott metal-insulator transition, 
and show that the sigma-model provides an adequate 
description of either one. 



1. The nonlinear sigma-model 



We return to Eqs . ( 2.20| ), ( 2.21 ), and use the repre- 
sentation, Eq. (3.24), of Q in terms of S and P. It is 
convenient to define a transformed ghost field A by 



A(x) =5- 1 (x)A(x)5(x) 



(3.34) 



In terms of these varia bles, the partition function for the 
replicated theory, Eq. ( 2.20 ), reads 



Z = I D[P] D[S] D[A] exp(A [S, P, A] + Tr In I[P]) . 

(3.35) 



Here I[P] is the invariant measure, Eq. ( 3.25 ), and the 
action A reads, 

A [S, P, A] = A dls [P] + Ant [SPS- 1 ] 

+ ^Tr In (G- 1 -iSAS- 1 ) 



dxti (A(x)P(x)) , 



(3.36) 
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with Adis and A lnt denned by Eqs. ( |2.23| ) and ( |2.24| ). 

Now we expand S, P, and A about their saddle point 
values. Alternatively, one could expand about the ex- 
act expectation values. However, the 'equation of state' 
for the latter is just the saddle point equation of state 
plus loop corrections, and for our purposes the difference 
is irrelevant. We therefore do not distinguish between 
saddle point values and expectation values, and denote 
both by (S), (P), and (A), respectively. As a result, we 
do not distinguish either between the saddle point Green 
function as defined in Eq. (2.36b), and the Green func- 
tion that contains the full expectation value of A as a self 
energy. From Eqs. ( |2.32 ) it follows that 



(S) = 1 <g> r , 



(3.37a) 



P)u = $u (to ® s ) \iuj ni - p 2 /2m + fi 



-*<A 12 ) 



N F (tti2 + «7l2) 



5 12 (to <8> s ) - N F (tti + 171) 



(3.37b) 



(Ai 



iu n 0l paa\ 
e \ r nn I 



and we write 
P=(P) 



AP 



A = (A) + AA 



(3.37c) 



(3.37d) 



Here we have defined two new matrix fields, tt and 7, for 
later reference. Their diagonal elements, tt\ and 171, are 
odd and even functions, respectively, of their Matasubara 
frequency arguments. For small frequencies, 



7Ti = sgnw ni 



(3.37e) 



while 71 is a real c onstan t that depends on the momen- 
tum cutoff in Eq. ( |3.37bD . 

Now we use the cyclic propert y of t he trace to write 
the Tr ln-term in the action, Eq. (3.36), in the form 



Tr ln(C?o 1 - SikS- 1 ) = Tr In (S^G^S - iA) 

= Tr In (G- 1 ) + Tr ln(l + G^S" 1 (d T S) 



-- GspS- 1 (V5) V + ^- G sp S 



2m~ sp 
-G S pi(AA) 



- 1 <V 2 S) 



(3.38) 



Notice that now the transformation matrix S always ap- 
pears in conjuction with some derivative, while the fluc- 
tuations AA are massive.Ea We proceed to expand in 



powers of the derivatives of S, and in powers of AA, in 
analogy to the procedure used in Ref. 18. The simplest 
term in this expansion is 



Tr (G sp S~ 1 (d T S)) = J dxtr (if2S(x)G sp (x 

x5- x (x)) 

J dxtT + ( o2 Q) 

Here we have defined a frequency matrix, 
^12 = (to <8> s ) S 12 aj ni , 



0) 



(3.39) 



(3.40a) 



and a new field, Q(x) — S(x) (tt + i7)6> _1 (x). Since 
the contribution of 712 to (P) 12 is, in the low-frequency 
limit, just an imaginary constant, it leads only to a con- 
stant contribution to Q, which does not contribute to 
the lowest order terms in the present expansion. We can 
therefore neglect it, and write 



Q(x) =5(x)^5- 1 (x) , 



(3.40b) 



with tt from Eq. (3.37b). Notice that Q is hermitian, and 
has the properties Q 2 = 1, and tr Q = 0. Q will turn out 
to be a convenient parametrization of the soft modes, in 
analogy to the unit vector (f> in </> 4 -theory, see Appendix 

No w co nsider the gradient terms in the expansion of 
Eq. ( 3.38 ). To this end, it is convenient to define the 
vector field, 



s(x)=5- 1 (x)(V5)(x) 



(3.41) 



with matrix valued components s l (x), (i — 1, . . . d). The 
expansion now proceeds in powers of s, or V<S. The term 
linear in s vanishes by symmetry. To quadratic order in 
s, both the next-to-last term under the Tr In (after a 
partial integration) and the square of the preceding term 
contribute. The result is a contribution to the action, 



7~T v 12 ^3 fa) / dx s i2 ( x ) s 2i ( x ) 



(3.42a) 



with 



?7i2,y(q) 



- [G sp (q,w ni ) + G sp (q,w„ 2 )] 
■ Qi Qj G sp (q, u ni )G sp (q, w„ 2 ) 



(3.42b) 



The approximation that neglects all fluctuations of P 
and A implies a factorization of four-point functions into 
products of Green fapctions, and hence Eq. (3.33b) ap- 
plies. We thus haveSj 

lim n t7 12 Vi2,ij{c0 = \ (1 - 7T17T2) % V , (3.43a) 

n 1 ,n 2 —>0 V z — ' Z 
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where 



(3.43b) 



with G± = G sp (uj n -> ±0). The structure of Eq. ( fr.43a| ) 
allows one to write the low-frequency limit of the term 
quadratic in s in terms of the matrix Q defined in Eq. 

(Hi), 

J2 [ dX S 12( X )4l( X ) h^2vi2,ij{^} 
12 q 

= ^-V [dx.tr (vQ(x)) 2 + 0(V 3 ,OV) . (3.44) 



Finally, we can use Eq. (3.33b) again to rewrite the cou- 
pling constant rj in a more familiar form, 



n 



-l 

2dmV 



X> 2 [G+(q)-G_(q)] 2 = 7rmc7o . (3.45) 



One recognizes <7q as the conductivity in the self- 
consistent Born approximation with the interaction taken 
into account in Hartree-Fock approximation. 

The remaining task is to deal with the fields P and A. 
Since the fluctuations of these fields are massive modes, 
in order to capture the leading effects of the soft modes 
in the system it suffices to integrate them out in any ap- 
proximation that preserves the Ward identities derived 
in Sec. Ill A. The simplest way to do so is to integrate 
them out in saddle-point approximation, i.e. to neglect 



AP and A A everywhere. From Eqs. ( 3.36| ), ( |3.37 ) we see 
that then the entire action can be expressed in terms of 
Q. Finally, we subtract from Q its saddle-point value 
(Q) = an d define a matrix field Q = Q — n. Replacing 
Q by Q everywhere in the action leads just to an uninter- 
esting constant contribution. This is obvious for all terms 
except for the r = channel in the spin-singlet interac- 
tion term. There a term linear in Q seems to remain, but 
inspection shows that it is proportional to tr Q = 0. We 
then obtain the following effectice action, 

AnloM = Co J dx tr ^VQ(x)j 

+ ^~r" / ^ ^ (^W) +Aint[Q] , 

(3.46a) 

in terms of the matrix field Q = Q — ir, where ir is the 
diagonal matrix defined in Eq. ( 3.37b| ), and Q is subject 
to the constraints, 



2 (x) = l®r , Q f = Q , trQ(x)=0 



(3.46b) 



We note that the hermiticity of Q is a consequence of our 
having dropped its anti-hermitian part, since it does not 
contribute to the leading terms in the long- wavelength, 
low-freq uency effective action. The effective action given 
by Eqs. (3^49) is the generalized nonlinear sigma-modcl 



that was first proposed by Finkel'steinEj as a model for 
interacting disordered electrons, and whose properties 
have been studied in considerable detail, in particular 
with respect to the metaL-insulator transition that is de- 
scribed by the model.ESEJ In the next subsection, we will 
discuss AnlctM, as well as its corrections, from a RG 
point of view, and show that it suffices for a correct de- 
scription of both the critical FP and the stable Fermi- 
liquid FP, as well as the leading corrections to scaling at 
the latter. 



2. The disordered Fermi-liquid fixed point 

In this subsection we discuss the RG properties of the 
model derived in the previous section and, in particular, 
introduce the concept of a disordered Fermi-liquid FP. 
Physically, this FP characterizes a system with a finite 
density of states at the Fermi surface and diffusive two- 
particle excitations. 

a. The fixed point action To proceed, it is convenient 
to parametrize the matrix field Q in a way analogous to 
the 7r-field parametrization of the O(N) vector model, 
see Appendix |b|. We thus write Q in a block matrix form 
analogous to that used in Eq. ( 3.23| ) as 



qq< 



<T 



(3.47) 



where the matrix q has elements q nm whose frequency 
labels are restricted to n > 0, m < 0. This rep resentation 
builds in the constraints given by Eq. ( 3.46b ). While the 
sigma-model can be entirely expressed in terms of Q, this 
is not the case for the corrections to it. We therefore also 
express S in terms of q. For what follows, we will need 
only the first terms of an expansion in powers of q. From 
Eqs. ( |3.40b| ) and ( fTI?] ) we obtain, 



S = 1 ® T 

= 1 <8> TQ 



AS 
1 



-q 



0(q 2 ) 



(3.48) 



Now we perform a momentum-shell RG procedure! 
For the rescaling part of this transformation, we need 
to assign scale dimensions to the s oft fie ld q, and to the 
massive fields AP, AA from Eq. (3.37d). Choosing the 
scale dimension of a length L to be [L\ = —1, we write 
in analogy to Eqs. (B6), 



[q(x)] = -(d-2 + r,') 



[AP(x)] = [AA(x)] = -(rf- 



V) 



(3.49a) 



(3.49b) 
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which defines the exponents r\ and ry'. The stable Fermi=j 
liquid FP of the theory is characterized by the choice,t3 



n 



rf = 



(3.50a) 



Physically, rf — corresponds to diffusive correlations of 
the q, and rj = 2 means that the correlations of the AP 
are of short range. This is indeed what one expects in 
a disordered Fermi liquid. In addition, we must specify 
the scale dimension of frequency or temperature, i.e. the 
dynamical scaling exponent z = [u>] = [T]. In order 
for the FP to be consistent with diffusion, that is with 
frequencies that scale like the squares of wave numbers, 
we must choose 

z = 2 . (3.50b) 



Now we expand the sigma-model action, Eq. ( |3.46a ), 
in powers of q. In a symbolic notation that leaves out 
everything not needed for power counting purposes, we 
write 

A N LaM = -^- J dx (Vq) 2 + H J dxtoq 2 

+TT J dxq 2 + 0(X7 2 q i ,ojq 4 ,Tq 3 ) , (3.51) 

with the bare coupling constants G ~ l/oo, and H ~ 
Np. r can stand for any of the three interaction coupling 
constants. Power counting shows that, with the above 
choices for the exponents, all of these coupling constants 
have vanishing scale dimensions with respect to our FP, 



[G] = [H] = [T] = 



(3.52) 



These terms therefore make up part of the FP action. 

Now consider first the corrections that arise within the 
sigm a-mo del. The leading ones have been indicated in 
Eq. ( |3.5l| ). We denote the corresponding coupling con- 
stants by Cv2 g 4, etc., with a subscript that identifies the 
structure of the respective contribution to the action. 
One finds 



(3.53a) 



With respect to the remaining term we note that it is of 
order q 3 . Any contribution to physical g-correlation func- 
tions therefore contain this term squared, and therefore 
the relevant scale dimension is 



ctv) 2 = 2 [c T?3 ] = -(d - 2) 



(3.53b) 



We see that all of these operators are irrelevant with re- 
spect to the Fermi-liquid FP for all dimensions d > 2, 
and that they become marginal in d = 2 and relevant in 
d < 2. All other terms that are contained in the sigma- 
model are more irrelevant than the ones considered. 

We next consi der c orrections to the sigma-model ac- 
tion. From Eq. (3.36|) we see that there are five classes 



of such terms: (1) Adis, (2) the term tr (AP), (3) those 
parts of the Tr ln-term that contain the massive fluctu- 
ation AA, (4) the parts of Amt that contain the massive 
fluctuation AP, and (5) co ntrib utions to (P ) that were 
neglected in writing Eqs. (3.39) and (3.401). Since all 



terms linear in the fluctuations vanish, class (1) just con- 
tributes 



A dls [AP} = - — I dx (tr AP(x)) 2 

n 



. JL f dx tr (AP(x)) 2 , 

Trel J 



(3.54) 



where we have scaled the relaxation times T\ and r re ; 
from Eqs. ( 2.23| ) with l/2nNp and 1/nNp, respectively. 
Power counting yields 



[n] 



[Trell 







(3.55) 



so Adis is part of the FP action. Class (2) also con- 
tributes a marginal term to the action, and so does the 
leading contribution of class (3), which has the structure 
Tr(GAA) 2 . Combining these two marginal terms, we 
have a contribution to the FP action, 

A AP = J dx tr (AA(x) AP(x)) 

+ i J dx dy tr (G(x - y) AA(y) G(y - x) AA(x)) , 

(3.56) 

apart from an uninteresting constant. The remaining 
terms of class (3) are irrelevant contributions that couple 
q and A A. The least irrelevant terms have the struc- 
tures Tr (AAwg), and Tr (AAV 2 q). For contributions 
to physical correlation functions we must again consider 
the squares of these structures, and the relevant scale 
dimensions are 



(3.57) 



These terms are more irrelevant than the corrections con- 
tained in the sigma-model for all d < 4. Moreover, due 
to the even integer value of the scale dimension, they lead 
only to analytic, and therefore uninteresting, corrections 
to scaling near the FP. The leading term of class (4) is 
one that couples q (from one of the S) and AP. The 
corresponding coupling constant has a scale dimension 









(cAAuiq) 2 




( C AA V 2 q) 



(cTqAp) = -2 



(3.58) 



Again, these operators lead only to analytic corrections 
to scaling. Cross correlations between these various 
terms also appear, and the corresponding operators have 
the same scale dimension, namely —2, as the ones shown. 
Finally, it is easily checked that the least irrelevant terms 
of class (5) are of the structure 0J 2 q 2 : and they also have 
a scale dimension 
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(3.59) 



Tv(T,u) = 1 v{Tb z 1 ub- < - d -V) 



(3.62a) 



The above RG arguments show that the theory con- 
tains a Fermi-liquid FP that is stable for all d > 2, 
and analogous to the low-temperature FP of an O(N) 
Heisenberg model. The FP action cons ists of the nonlin- 
ear sigma-model action AnlctM, Eq. ( 0.46a ), exp anded 
to se cond order in q, and of Adis and A\p, Eqs. (3.54) 
and (3.56). We have further shown that the leading non- 
analytic corrections to scaling near the Fermi-liquid FP 
are contained in the nonlinear sigma-model. An effec- 
tive action that describes both the FP and the leading 
corrections to scaling is therefore 



A ef f = A 



NLcrM 



q] + A dls [AP] + „4 A p[AA, AP] 



(3.60) 



Notice that the soft and massive modes do not couple in 
this effective action. Furthermore, the massive fields AP 
and AA appear only quadratically. If desired, then the 
latter can be integrated out to obtain an action entirely 
in terms of the physical fields q and AP. 

b. Scaling behavior of observables We now discuss 
the physical meaning of the corrections to scaling in- 
duced by the irrelevant operators that we have identi- 
fied above. Let us denote by the generic name u any 
of the least irrelevant operators whose scale dimension is 
[u] = — (d — 2), and let us discuss various observables, 
viz. the conductivity a, the specific heat coefficient jy, 
the single-particle density of states N, and the spin sus- 
ceptibility Xs ■ Which of the various operators with scale 
dimension — (d — 2) is the relevant one depends on the 
quantity under consideration. 

Let us first consider the dynamical conductivity, a{ui). 
Its b are v alue is proportional to 1/G, and according to 
Eq. ( 3.52 ) its scale dimension is zero. We therefore have 
the scaling law, 



a(w,u) =o-(Lub z ,ub-( d ~^) 



(3.61a) 



where b is an arbitrary RG sc ale fac tor. By putting b = 
1/oW 2 , and using z = 2, Eq. ( 3.501 ), as well as the fact 
that o~(l,x) is an analytic function of x, we find that 
the conductivity has a singularity at zero frequency, or a 
long-time tail, of the form 



o~(oj) — const. 



,(d-2)/2 



(3.61b) 



This nonanalyticity is well known from perturbation the- 
ory for both noninteractingC-3 and interactingcJ electrons. 
This shows that i n this case either one of the coupling 
constants in Eqs. ( 3.53a ) and ( 3.53b ) can play the role of 
u. The present analysis proves that the w^ -2 )/ 2 is the 
exact leading nonanalytic behavior. 

The specific heat coefficient, ■yv — Cy /^niffri? ro P or " 
tional to the quasiparticle density of states ilffla whose 
scale dimension vanishes according to Eq. ( 3.52 ). We 
thus have a scaling law 



which leads to a low-temperature behavior 

7v (T) = const. + T^-V' 2 . (3.62b) 

From perturbation theory it is knownH that 7v shows 
this behavior only for interacting electrons, while for non- 
interacting systems the prefactor of the nonanalyticity 
vanishes. In this case, therefore, u must be equated with 
i c Tq 3 ) 2 - This can not be seen by our simple counting 
arguments. 

The single-particle density of states, N, is proportional 
to the expectation value of Q, and to study the leading 
correction to the finite FP value of N it suffices to replace 
Q by Q. Then we have, in symbolic notation, N ~ 1 + 
(iQ ) + • • ■ = 1 + AN. The scale dimension of AN is 
[AN] = 2[q] = d - 2. We find the scaling law 

AN(u}) = b^ d - 2) AN(uub z ) , (3.63a) 

which leads to the so-called Coulomb anomaly^! 

N{u) = const. W d ~ 2)/2 , (3.63b) 

Again, this behavior is known to occur only in the pres- 
ence of electron-electron interactions. 

Finally, we consider the wave vector dependent spin 
susceptibility, Xs{<\)- Xs is given by a Q-Q correlation 
function, and the leading correction to the finite Fermi- 
liquid value is obtained by replacing both of the Q by q. 
Then we have a term of the structure \s ~ T J dx (q^q), 
with scale dimension [% s ] = 0. The relevant scaling law 
is 



Xs(q,u) = Xs (qb,ub-( d - 2 )) 



(3.64a) 



which leads to a nonanalytic dependence on the wave 
number, 



X s (q) = const. + |q 



(d-2) 



(3.64b) 



This behavior is also known from perturbation theory, 
and holds only for interacting electrons. It has recently 
been shown to have intetesting consequences for the the- 
ory of ferromagnetism.113 

To summarize, we see from the above argupijcnts 
that all of the so-called weak-localization effects,E3 i.e. 
nonanalytic dependencies of various observables on fre- 
quency, temperature, or wave number in disordered elec- 
tron systems that are well known from perturbation the- 
ory, emerge naturally in the present context as the lead- 
ing corrections to scaling near the Fermi-liquid FP of a 
general field theory for disordered interacting electrons. 
Apart from providing an aesthetic, unifying, and very 
simple explanation for these effects, our arguments also 
prove that they do indeed constitute the leading non- 
analytic behavior, a conclusion that cannot be drawn 
from perturbation theory alone. We emphasize the ease 
with which these results are obtained within the present 
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framework. This reflects the judicious choice of our start- 
ing point for the RG analysis. The transformation from 
fermionic variables to composite ones, together with the 
symmetry analysis that identifies the soft modes, pro- 
vides for a formulation of the theory that allows for rather 
nontrivial results to be obtained from a simple RG anal- 
ysis at the Gaussian level. Similar benefits can be ob- 
tained from analogous treatments of other systems, see 
Appendix |b|. 

Finally, it is interesting to note that similar nonana- 
lyticities occur in very different systems, for example, in 
classical fluids. In that context they are known as long- 
time tail effects, and were first discussed theoretically by 
using many body perturbation theory and mode coupling 
theory.EJ Later, they were examined using RG ideas, and 
they were shown to be related to corrections to the scal- 
ing behavior near a hydrodynamic FPE3 



3. The critical fixed point 

In addition to the stable Fermi-liquid FP discussed 
above, the ge nerali zed matrix nonlinear sigma-model 
given by Eq. ( 3.46a ) also possesses another FP that in 
general has the properties of a critical FP. It characterizes 
a zero-temperature metal-insulator transition in dimen- 
sions greater than two. For a review of this FP, and of the 
various universality classes for the metal-insulator tran- 
sition, we refer the reader elsewhere.t£l Here we briefly 
comment on the differences between the Fermi-liquid FP 
and the critical one, and on the justification for using the 
sigma-model to describe the critical behavior. 

First of all, at the critical point the spontaneously bro- 
ken symmetry is restored, so the distinction between pos- 
itive and negative frequencies vaaishes as the magnitude 
of the frequencies goes to zero.Eil Physically this corre- 
sponds to a vanishing single-particle of density of states 
at the Fermi level, which plays the role of an order pa- 
rameter. This characterizes the metal-insulator transi- 
tion as being of Ander son-M ott type. For example, the 
left-hand side of Eq. (3.17a) vanishes at the transition 



as ni,ri2 — > 0. Therefore, matrix elements Q ni n 2 with 
n\ri2 > and n\ni < 0, respectively, should both scale 
in the same way. This implies that at the critical FP, rj' 
in Eq. (3.49a) is given by 



n 



2-d. 



(3.65) 



so that q is dimensionless. The physical critical exponent, 
on the other hand, is given by 



V 



(3.66) 



To see this, we note that Q is proportional to (P), which 
in turn is proportional to the density of states which van- 
ishes at the transition with a critical exponent /3H Here 
v is the critical exponent that describes the divergence 
of the localization or correlation length £. It occurs in 



the scaling equality, Eq. ( 3.66| ) , since r\ characterizes mo- 
mentum singularities at criticality. 

We add some comments on the scaling of AP fluctu- 
ations. At the critical point, the AP fluctuations are 
actually of long range, and the sigma-model description 
breaks down. To avoid this problem, one must work at 
sufficiently small momenta p (or frequencies) , and at suf- 
ficiently large distances t from the critical point, so that 
the q fluctuations still dominate over the AP fluctua- 
tions. Since the AP correlation function is related to the 
longitudinal susceptibility, it diverges as i~ 7 while the q 
correlations diverge as p~ 2+r >£3 The relevant inequality 
is therefore 



(3.67) 



where we have used the scaling equality v = 7/(2 — rf). 
This can also be seen more explici tly b y considering the 
AP-AP propagator. From Eq. ( |3.54| ), and the corre- 
sponding term that is proportional to (V(AP)) 2 one sees 
that the AP- mass mp has a scale dimension of 1 (see 
also Appendix B3). The AP fluctuations can therefore 
be neglected for momenta that ar e sma ll compared to 
in p r~j ~ t u , which is again Eq. (3.67). In this regime 
the AP fluctuations are effectively of short range, just as 
they are at the Fermi-liquid FP, and the sigma-model is 
valid. Even though the range of validity of the sigma- 
model goes to zero as the critical point is approached, 
it still can be used to extract the critical behavior. The 
salient point is that the RG, and the scaling behavior de- 
rived from it, allow s one to extrapolate from the regime 
given by Eq. (3.67) to the critical region. However, no 
information can be obtained in this way about the sym- 
metric phase, i.e. the insulator. 

Apart from these theories that work near d = 2, 
the matrix nonlinear sigp.a-model has also been stud- 
ied in high dimensions .cJ This work has established 
rf+ = 6 as an upper critical dimension, above which the 
Anderson-Mott transition is correctly described by a sim- 
ple Landau-type theory. This treatment of the problem 
stresses that the metal-insulator transition problem is, 
somewhat counterintuitively, conceptually simpler in the 
presence of electron-electron interactions than in their 
absence, since the interacting problem possesses a sim- 
ple critical order parameter, viz. the density of states 
at the Fermi level. Also, in these papers the presence of 
random-field like terms in the renormalized action was 
discovered. These terms are responsible for the upper 
critical dimension being 6 (rather than 4 as in the O(N) 
vector model), and they have led to the suggestion that 
the Anderson-Mott transition may have features remini- 
scient of a glass transition.EJ 



4- Ferromagnetic transition in the metallic phase 



It has been known for some time that the matrix non- 



linear sigma-model, Eq. (3.46a), contains another critical 
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FP that is not related to a metal-insulator transition.ES 
While originally the nature of the corresponding phase 
transition was not clear, we have recently shown that 
it is the zero-tcmpcrature transition from a paramag- 
net to a ferromagnet ip-ia disordered metal in dimen- 
sions greater than two.E3 In this reference, a Landau- 
Ginzburg- Wilson (LGW) functional for spin density, or 
order-parameter, fluctuations has been derived. To ob- 
tain this functional, we integrated out all excitations or 
modes other than the order parameter, including all of 
the soft mo des related to weak localization effects, see 
Sec. Ill B 2 above. In the effective, or order-parameter, 



field theory these extra soft modes led to nonanalyti- 
cies in the bare LGW functional. Using rcnormaliza- 
tion group methods, we were able to exactly determine 
the critical behavior at the magnetic phase transition.ta 
The earlier approach of Ref. |55|, on the other hand, had 
focused entirely on the behavior of the diffusion modes 
across the magnetic phase transition, while the order pa- 
rameter fluctuations had effectively been integrated out. 
In some respects, these two approaches to the problem 
were therefore complementary to one another. 

Even though both of these theories led to the same 
results, their physical underpinnings are very different, 
and neither one of them constitutes what one might con- 
sider the physically most obvious approach. Physically, 
the most sensible procedure would be to treat all of the 
soft modes at this quantum phase transition on the same 
footing. These soft modes include the diffusive modes 
discussed above, i.e. the Q nin2 with n\ri2 < 0, as well as 
the fluctuations of the magnetic order parameter, that is, 
the spin density or magnetization. In terms of Q-fields, 
the local spi n dens ity at frequency u n , n l s (x, u) n ), is given 
by (cf. Eq. fl2.12b|) 



m r— 0,3 



(-l)"C + nW > (3-68) 



and the macroscopic magnetization order parameter, M, 
has components, 



M' 



= (n>„ = 0,x)) 

= ^££(-d t 



n ( \ 



(3.69) 



m r=0,3 



From these equations we see that, as the magnetic tran- 
sition is approached, the c lassification of soft modes that 
was given in Section III A above, breaks down. New soft 
modes related to the criticality of Q m m occur. Alterna- 
tively, the extra soft mode that appears near the ferro- 
magnetic transition can be related to a new zero eigen- 
value in th e Gau ssian eigenva lue problem. The matrix 
M in Eqs. ( 2.37a| ) and ( |2.39a ) is proportional to, 



* M 

711112,713114 



(P) 



5 n 1 n 3 5 n 2 n 4 2? ra 1 n 2 



(P) 



-2TT t 5. 



(3.70) 



with i = 1,2,3 and we have considered the term diagonal 
in all of the replica labels since these terms contain the 



interactions that lead to magnetism. The corresponding 
eigenvalue problem is, 



n 

= X fn in2 (P) , 



(3.71) 



with / the eigenfunction and A the eigenvalue. Setting 
ni = ri2 and integrating over n\ leads to an equation for 
A, 



1 = 2r t r 



P»n(p = 0) 

l-AD„ n (p = 0) 



(3.72) 



Criticality is characterized by a new zero eigenvalue, A = 
0, and we have put p = 0, si nce t his is where the first 
zero eigenvalue occurs. Eq. (3.72) reduces to a Stoner 



criterion for the occurrence of ferromagnetism, modified 
by disorder. Expanding in powers of A leads to 



A 



-t 



(3.73) 



with t the dime nsionless distance from the critical point. 
From Eq. (3.71), we see that the critical eigenfunction is 



/; m (p) = 2> nn 2r t vTn>» - 0,p) 



(3.74) 



The above discussion makes it clear that a physi- 
cally satisfactory theory of the ferromagnetic, or any 
other, phase transition in a disordered metal should 
take both the soft modes that were discussed in the 
preceding subsections, and the additional, critical, soft 
modes into account on the same footing. Techni- 
cally one needs to extract the critical part from the 
P-fluctuations. The details of such a theory of the 
paramagnetic-to-ferxomagnetic phase transition will be 
discussed elsewhereH] Here we just mention why the 
sigma-model, in Refs. g5|, manages to describe the correct 
critical behavior, even though it neglects the soft modes 
related to the magnetization. The salie nt po int is again 
that there is a regime, described by Eq. ( |3.67 ), where the 
diffusive modes that are contained in the sigma-model 
dominate over the magnetization fluctuations that are 
not. As in the case of the metal-insulator transition, the 
RG allows one to extract the critical behavior from an 
analysis of that region. 



IV. THE CLEAN LIMIT 

In this section we discuss the clean limit of our field 
t heor y. For t he re asons that were pointed out after Eqs. 
( |2.13 ) in Sec. HA, this treatment of clean Fermi systems 
will not be as complete as our theory for the disordered 
case. For notational simplicity we will also suppress the 
particle-particle degrees of freedom, but they can be eas- 
ily restored. Even with our restriction to two Fermi- 
liquid parameters (in the absence of the Cooper channel) 
instead of infinitely many, we will be able to study funda- 
mental structural properties of the clean limit that will 
have to be included in any more complete theory. 
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A. Effective Q-field theory 



Let us return to the action, Eq. ( 2.21 ), and perform 
the clean limit, ri,r re ; — > oo. As pointed out earlier, 
the theory has been purposely set up so that this limit 
can be taken without diffic ulty. Adis then vanishes. The 
saddle point solution, Sec. II C , reduces to the ordinary 
Hartree-Fock approximation, and the Gaussian approx- 
imation reduces to RPA. In particular, the saddle point 
action contains an exact description of the noninteracting 
electron gas. This will be important in what follows. 

As was pointed out at the end of Section |fi|, the fluc- 
tuations of A are not massive in the clean case, since 



the function tp, Eq. ( 2.38b ), is singular in the long- 
wavelength, low-frequency limit, see Eq. ( [2.43 ). In- 
deed, Eqs. ( |2.37 ) and ( |2.39 ) show that the Gaussian A- 
propagator is just minus the noninteracting part of the 
Q-propagator. On the general principle that one should 
keep all of the soft modes on the same footing, this seems 
to suggest keeping the A field at the same level as the Q- 
field. Dropping the constant sad dle-p oint contribution, 
and expanding the Tr In in Eq. (2.21), the action then 
reads 



A[Q,A] = A C Q n [Q,A] + AA[Q,A] 



(4.1a) 



where A C Q ean = Ac — Adis with Ac and Adis horn Eq. 



(2.36a), and 



AA[Q,A] 



— 2 M ~ 1 i M 



E 



E E 



M y( M - 2 )/ 2 

M=3 ki,...ki( ni...fm 

x ^ki + ...+k M ,0 X^.„ M ( kl > ■ ■ ■ > k M-l) 
x VVn 2 ( k l)'''^nin 1 ( k M) 

xtr [(((5A) ni „ 2 (ki) - (<5Q) ni „ 2 (ki)) • • • 

(kM))] ■ 



(4.1b) 



Here tp is given by Eq. ( 2.38b[ ), and 



V 



xG sp (p + ki, LOn 2 ) ■ ■ ■ G sp (p + kjlf-1, w„ M ) 



(4.1c) 



In this formulation of the theory the vertices are given 
by the product of the M-point corr elati on function x^ M ^ 
and the M f actors of in Eq. ( 4.11 ). An inspection 
of Eq. (4.1c) shows that the generic behavior of > 
considered as a function of some generic wave number 
k, is characterized by a divergence l/fc M_1 for small k. 
This is the most divergent behavior ^( M ) can display, 
and it is realized unless either all of the frequencies ui ni 
through uj nM are positive, or all of them are negative, 
in which case x scales like a constant. The behav- 
ior of the product of the p^ 1 , on the other hand, de- 
pends on the detailed distribution of the frequency la- 
bels. If riirii+i > 0, then Pn 1 n +l i s a number, while for 



nifii+x < it goes like k for small k. If only two of the 
frequency pairs have elements with opposite signs (this 
is the smallest nonzero number possible), then the com- 
plete vertex scales like l/fc M_3 . For M > 3, the vertices 
are thus not finite in the limit of long wavelengths and 
small frequencies. In other words, the Q-A field theory 
given by Eqs. (4.1) is not local. The reason for this is as 
follows. As was mentioned in Sec. ||, in the clean limit 
the one-particle excitations are soft, as is manifested by 
the massless single-particle Green function. These soft 
modes have been integrated out when we integrated out 
the fermions, and this is what leads to the nonlocalities 
in the effective field theory. The massless single-particle 
excitations are also indirectly responsible for the softness 
of the A fluctuations, and for the need to keep infinitely 
many Landau parameters. All of these difficulties thus 
have the same underlying source. 

Given that the matrix field theory is nonlocal in any 
case, we can proceed and integrate out A, as it will turn 
out that this does not lead to further undesirable prop- 
erties of the theory. The integrating out of A can be 
done exactly in the sense of a prescription for doing per- 
turbation theory for the resulting Q-field theory. To see 
this, we recall that the A propagator is minus the non- 
interacting part of the Q-propagator. As a result, the 
integration over A just cancels the noninteracting parts 
of any internal Q-propagators in a loop expansion for any 
Q-correlation function. This is easy to see for the first 
few terms in the loop expansion, and it can be proved by 
induction to be true order by order in perturbation the- 
ory. As an illustrative example, the cancellation scheme 
is demonstrated diagrammatically in Figs. || and|^ for the 
two-point vertex function, and the two-point propaga- 
tor, respectively, to one-loop order. Notice that, in order 
to avoid double counting, one effect of the A-field must 
be the cancellation of the noninteracting contributions, 
since the saddle-point contribution to the action already 
contains a complete description of the noninteracing elec- 
tron system. What is remarkable is that it does nothing 
else. We then obtain the following effective action en- 
tirely in terms of Q-fields, 



A[Q] = A a [Q] + Am [Q] + AA[Q] , 



(4.2a) 



where 



k nm r 7 i 

x;W)„ m (-k) , (4.2b) 
is the noninteracting part of the Gaussian action, and 

A int [Q] = -8]T ]T r«£T ]T r (SQ) nir Jk) 

tO r=0,3 k ni,ri2,n3,ri4 

x<5 ni _ n2 ,„3_„4(<5Q) n3n4 (-k) , (4.2c) 

is the interacting one. The (i = 0, 1,2,3) were de- 
fined after Eq. ( fP9b| ). The non-Gaussian part of the 
action reads, 
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+ (2-loop terms) 





+ (2-loop terms) 

FIG. 2. Perturbation theory to one-loop order for the 
2-point Q-vertex function. The small circle denotes the Gaus- 
sian vertex, the solid and dashed lines are Q-propagators and 
A-propagators, respectively, and the solid lines with a vertical 
bar denote the interacting part of the Q-propagator. 



AA[Q] 



g 2 M - J {-i) M 



E E 

M—3 ki,...k Jf n 1 ...n M 



X(5 kl + ...+k M ,0 Xnf)..n M 



XV ni n 2 ( k O---Vnini( k M) 

xtr [(5Q) (ki).-.(tfQ) nMBl (k M )] • (4.2d) 

This effective action needs to be supplemented by the 
following rules for doing perturbation theory: 

rule 1: For calculating Q-propagators, all internal 
propagators must be taken as the interacting part of 
the Gaussian propagator, i.e. as the second term on the 
right-hand side of Eq. (2.39a). 

rule 2: For calculating Q-vcrtcx functions, rule 1 also 
applies. In addition, one needs to consider all reducible 
diagrams (which normally do not contribute to the ver- 
tices), with all reducible propagators replaced by minus 
the noninteracting Gaussian Q-propagator, i.e. minus 
the first term on the right-hand side of Eq. ( 2.39a| ). 

Explicit calculations using these rules readily show 
that the integrals that correspond to the diagrams in 
a loop expansion are identical to integrals that one en- 
counters in standard many-body perturbation theory for 
the same quantity. This allows for a one-to-one corre- 
spondence between many-body diagrams and the loop 
expansion based on the present field theory. Neverthe- 
less, even at a calculational level the present formulation 
provides advantages compared to standard perturbation 
theory. For instance, due to the above rules, the loop 
expansion is equivalent to an expansion in powers of the 
screened Coulomb interaction, with the zeroth order, i.e. 
the Gaussian theory, reproducing RPA. The loop expan- 
sion therefore allows for a systematic improvement over 
RPA. Perhaps more importantly, our field-theoretical 
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FIG. 3. Perturbation theory to one-loop order for the 
2-point Q-propagator. The notation is the same as in Fig. 
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formulation allows for an application of the renormal- 
ization group to draw structural conclusions about the 
theory in a nalogy to those discussed for the disordered 
case in Sec. IIIB 2. This we will discuss next. 



B. The Fermi— liquid fixed point 

We now are looking for a stable RG fixed point that 
describes the clean Fermi liqui d, in an alogy to the disor- 
dered Fermi-liquid FP of Sec. [II B 2. For this purpose 
it is again convenient to split the matrix Q into blocks in 
frequency space, 



Qr 



P 



,,,„ W inn • 
otherwise 



(4.3) 



As in Sec. Ill B 2 , we define the scale dimension of a 
length £ to be [L] = — 1, and we define exponents r\ 
and r/' by writing 



[?(*) 



[AP(x) 



(d - 1 + ?/) 



1 



(d - 1 + tj) 



(4.4a) 



(4.4b) 



Here AP = P — (P), and as in Sec. |irj we do not dis- 
tinguish between AP and SP. The FP action has the 
properties one expects from a Fermi liquid if we choose 



V 



1 



?/ = 



and a dynamical exponent 

2 = 1 



(4.5a) 



(4.5b) 



Power counting shows that with these choices, Ao is di- 
mensionless, and hence part of the FP action. So is the 
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part of Ai n t that is quadratic in q. The parts ofAmt that 
couple q with AP, and AP with itself, are irrelevant with 
scale dimensions of —1/2 and —1, respectively. 

Now consider the term of Q( Q M ) in the non-Gaussian 
part of the action, AA, Eq. ( |4.2d| ). Let N of the Q- 
fields be P's, and M — N be q's. Denote the coupling 
constant for that term by un.m-n- Taking into account 
the properties of x^ M ' an d <P discussed after Eqs. (4.1), 
we obtain 



[um.o] =-(d+l) 



(M - 2) 



(4.6a) 



for N = M, and 



[un,m-n\ 



i (M — A' — 2) — (cf — 1) 2~~ ^ 



(4.6b) 



for N 7^ M. Here we have made use of the following 
observation. Because of the two rules in the preceding 
subsection, all contractions that are performed in AA to 
obtain renormalizations of the Gaussian action result in 
interacting propagators. Since the interaction carries a 
factor of temperature, see Eq. (4.2c), this results effec- 
tively in (M — 2) /2 factors of T that need to be taken 
into accoun t in a ddition to the terms that show explic- 
itly in Eq. (|Dd|). This factor of t( m ~ 2 )/ 2 corresponds 
to the M — 2 internal Q-fields that get contracted in a 
renormalization of the Gaussian term. We note that this 
feature is automatically built into the theory only after 
the A- field has been integrated out. Doing so is there- 
fore advantage ous. Since M — N is necessarily even, we 
see from Eqs. fl4.6| ) that all of the non-Gaussian terms 
are RG irrelevant with respect to the Fermi-liquid FP, 
provided that d > 1. For d = 1 there is an infinite set 
of marginal operators, which signals the instability of the 
Fermi-liquicLground state against the formation of a Lut- 
tinger liquidEj For d < 1 the Fermi liquid is unstable due 
to infinitely many relevant operators, as expected. 

We mention that there is no consensus on whether 
or not interacting Fermi systems in d > li-ye neces- 
sarily Fermi liquids. In particular, AndersonEH has pro- 
posed that there exists a stable FP that corresponds to 
a Luttinger-type liquid, at least in certain 2-d models. 
While we do not find such a FP, we stress that the above 
considerations do not constitute a proof that none exists. 
All we have shown is that the assumption of a Fermi- 
liquid FP in d > 1 leads self-consistently to the conclu- 
sion that this FP is stable. ThiSj-is-in agreement with 
a variety ofrOther RG arguments, BO and with explicit 
calculations xB However, we do not know what the basin 
of attraction for the Fermi-liquid FP is, and we cannot 
exclude the existence of other fixed points. 

Obviously, the discussion of the corrections to scal- 
ing and their consequences for the behavior of thermo- 
dynamic quantities that was given for the disordered case 
in Sec. IIIB 2 can be carried over. The only difference is 



which we again denote by u, is now [u] — — (d— 1). As an 
explicit example, let us consider the spin su sceptib ility 
Xs- It obeys a scaling law in analogy to Eq. ( |3.64a ), 

Xs (q,«) = X s (q6,^ (d - 1) ) , (4.7a) 

which leads to a nonanalytic dependence on the wave 
number, 



X s (q) = const. + |q 



d-l 



(4.7b) 



that the scale dimension of the least irrelevant operators, 



This is the leading wave number dependence of \s for 
1 < d < 3, and the leading nonanalytic one in all di- 
mensions. In d = 3 there is a logarithmic correction to 
scaling, which our power counting arguments are not sen- 
sitive to, and the behavior is q 2 In |q|. This behavior has 
recently beeit-pbtained by means of explicit perturbative 
calculations,^ and its implications for the paramagnet- 
to-ferromagnet transition at zero temperature have been 
discussed.E3 The nonanalytic behavior of other thermo- 
dynamic quantities, and of the quasiparticle inelastic life 
time, can be understood by means of analogous argu- 
ments. 



V. CONCLUSION 

In this paper we have given a general method to study 
the long- wavelength, low-frequency behavior of many- 
fcrmion systems, both with and without quenched dis- 
order. The crucial ideas are to first identify the slow 
modes of the system by using a symmetry analysis, then 
to separate these soft modes from the massive ones, and, 
finally, to use renormalization group ideas to eliminate 
the degrees of freedom that are irrelevant in the long- 
wavelength limit. 

Using these ideas we have accomplished a number of 
things. We first established that in a disordered system, 
a stable disordered Fermi-liquid FP is possible for d > 2. 
We showed that the so-called, weak localization effects in 
itinerant electronic systemsE3 are, or can be interpreted 
as, corrections to scaling near this FP. This derivation 
not only reproduces known perturbative results, but also 
establishes that their functional form is asymptotically 
exact in the long- wavelength limit. In this respect our 
achievement is analogous to that of Rcf. |50| for long-time 
tails in classical fluids. In Section HI B| we have given 
a technically satisfactory derivation of the generalized 
nonlinear sigma-model that has been used to describe 
metal-insulator transitions near two dimensions. Finally, 
we have indicated how the theory must be modified near 
other quantum phase transitions where additional soft 
modes, namely the critical modes, appear. 

For fermion systems without disorder, analogous re- 
sults have been obtained. First, in agreement with oth- 
ers, we find that a clean Fermi-liquid FP exists for d > 1. 
Corrections to scaling near this FP show that clean elec- 
tronic systems have nonanalyticities in various correla- 
tion functions that are analogous to weak-localization 
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effects in disordered systems. This is again in agree- 
ment with results obtained on the basis of either Fermi- 
liquid theoryja or many-body perturbation theory. Our 
treatment provides a unified description that reveals a 
deep connection between the behavior of fermionic sys- 
tems with and without disorder, respectively. In Section 
1VA we also gave a novel perturbation theory method 
for clean fermion systems. 

There are still a number of things to be done. For 
disordered systems, we recently suggested an order 
rameter description of the metal-insulator transition 
Our theory was based on the nonlinear sigma-modcl ap- 
proach, applied to high dimensions (near d = 6), where 
its validity is not clear. This approach indicated that the 
metal-insulator transition has features in common with 
the transition in a classical random-field Ising model. 
These ideas need to be reexamined using a more general 
theory, since the sigma-model approach is asymptotically 
exact only for d < 4. Instanton solutions of the general 
field theory should also lead to insights concerning rare 
events, like local moment formation, and the effects of lo- 
cal moments on the metal-insulator transition. Finally, 
as already noted, additional quantum phase transitions 
from the Fermi-liquid state need to be investigated. 

For clean systems, the most important thing to do is 
to include the effects of other Fermi-liquid parameters, 
as was discussed in Section |n]. Once this is done, quan- 
tum phase transitions in clean itinerant systems can be 
properly studied. It will then also be possible to study 
the crossover from the clean to the disordered Fermi liq- 
uid fixed point in detail. For non-interaction systems, 
this latter poiat has been addressed by Muzykantskii and 
KhmelnitskiiH 

For both clean and disordered systems, it is interest- 
ing to ask how a Fermi-liquid FP can be avoided. This 
question has been of great interest recently in connection 
with high-T c superconductors, and other systems that 
have 'strange metal' phases. For clean systems this has 
recently been reviewed in Ref. |6^. For disordered sys- 
tems, the situation is less clear. Various itp>es of Kondo 
lattice mechanisms have been proposedJ£3 In these ap- 
proaches it is unclear how these effects modify the long- 
wavelength transport properties. Another possibility is 
to consider systems with a vanishing single-particle den- 
sity of states at the Fermi surface. If this occurs, then 
the soft modes discussed in Sec. [n] are not as singular, 
and a two-dimensional disordered metal phase becomes 
possible. Further, it suggests the possibility of an exotic 
metal-insulator transition in two dimensions. Clearly, 
these problems require much more work. 
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APPENDIX A: PROJECTION ONTO THE 
DENSITY 

Here we demonstrate, for the spin-singlet particle-hole 
interaction term, the proj ection onto density modes that 
was used in Sees. II A and II B 1 . An analogous procedure 
can be applied to the other interaction channels, and to 
the disorder part of the action. 

It is most convenient to go back to a Hamiltonian 
description of the system. The part of the interaction 
Hamiltonian that corresponds to Eq. (2.10b) reads 



ffil = J EE' r 2^) A/ k (q)A /p (q) 



k,p q 



where 



/k(q) = E c Ik c ^ 



(Ala) 



(Alb) 



is the phase space density operator in terms of electron 
creation and annihilation operators c' and c, respectively, 
and A(c t c) = cte — (c'c). In the space of products of 
fermion operators, we define the Kubo product,E-3 



(A\B) = dr (AB(t)AA* 
Jo 



(A2) 



with A and B operators, and r the imaginary time vari- 
able. In terms of this scalar product in operator space, 
the desired projector reads 



V 



in(q)) 



1 



K»(q)l 



(A3) 



Here n n (q) = ^ k /k(q) is the electron number density 
operator, and g(q) — (n n (q)\n n (q)) is the static density 
susceptibility or wave vector dependent compressibility. 
Using V twice, it is now easy to project onto the density 
in Eq. ( |Ala| ). We obtain 



H 



2^ 
q 



r (s) (q) n „(q) n n (- 



(A4a) 



where, 



r (s) (q) = -4rE r S(9)(/ k (q)K(q)) 

y ^ ' k,p 

x(n„(q)|/ p (q)) . 



(A4b) 



The phase space Kubo function, ,gk P (q) = (/k(q)|/ P (q)), 
for clean, free electrons is proportional to <5k P <5(k 2 — k 2 



so in this case Eq. ( A4b| ) results in pinning k and p 
to the Fermi surface. In a disordered system, <?k P (q) 
has a width given by the inverse elastic mean-free path, 
and hence r( s )(q) is a weighted average over a region in 
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the vicinity of the Fermi surface, as mentioned in Sec. 
II A . Finally, for physics that is controlled by soft modes 
and long-wavelength processes we can replace r^^(q) by 
r( s ) = rW(q = 0). Switching back to a field -theor etic 
representation of the fermions, we obtain Eq. (2.13a). 



APPENDIX B: O(N) SYMMETRIC 4 -THEORY 

In this appendix we perform an analysis of 4 -theory 
that is analogous to that of the matrix field theory in 
Sec. IIIB. Much of this material can be found in the 
literature, e.g. in Zinn-Justin's book, Ref. [35]. However, 
we find it useful to include it here for pedagogical rea- 
sons. The remarkable analogy between field theories for 
electrons apd classical spin models that was first noted 
by WegnertfJ substantially simplifies the understanding 
of the former in terms of the technically much simpler 
structure of the latter. 



1. Origin of the nonlinear sigma— model 

Let us consider an O(N) symmetric 4 -theory with a 
magnetic field h in the 1-direction. The action 



dx 



r(0(x)) 2 +c(V0(x)) 2 



+u j dx (^(x) -<Mx)J -hi dx(f>i(x) , (Bla) 

determines the partition function 

Z[h] = J D[$] e- s[ t ] . (Bib) 

In the low-temperature phase, where the O(N) symme- 
try is spontaneously broken, it is convenient to decom- 
pose the vector field (j) into its modulus p and a unit 
vector field 6, 



(x) = p(x) 0(x) , ^(x) = l 



(B2) 



(f> parametrizes the unit sphere, and thus provides a rep- 
resentation of the homogeneous space 0(N)/0(N — 1). 
In terms of p and <p the action reads, 



S \f, 4>\ = J rf x 
+ rp 2 (x) 



c^p 2 (x) (V0(x)) 2 +c< 2) (Vp(x)) 2 
-u J dxp i (x) — h J dxp(x)<j>i (x) , 



(B3a) 



and switching from the functional integration variables < 
to (p, (j)) leads to a Jacobian or invariant measure 



(B3b) 



In Eq. (B3a) the bare values of c^> and are equal, 
and equal to c. Notice that the field <f> appears only in 
conjunction with two gradient operators. <f> represents 
the N — 1 soft Goldstone modes of the problem, while p 
represents the massive modes. Now we parametrize <p, 



(x) = (a(x),f?(x)) , 



where 



CX(X) = v/l-7f2( X ) 



(B4a) 



(B4b) 



We split off the expectation value of the massive p-ficld, 
p(x) = M + Ap(x), with M = (p(x)), and expand in 
powers of 7? and Ap. Rescaling the coupling constants 
with appropriate powers of M, the action can be written 



S[p, 7?] = S NLaM [A + AS[p, 7f] 



(B5a) 



Here 
Snlom[it] 



dx 



(V7?(x)) 2 + (Va(x)) 2 



(B5b) 



—h J dx<j(x) 

is the action of the 0(N) nonlinear sigma-model, and 
A%,7r] = r J dx (Ap(x)) 2 + C W J dx (VAp(x)) 2 
+u 3 J dx (Ap(x) f + Ui J dx (Ap(x)) 4 



O Apa, Ap 



(B5c) 



If we neglect all fluctuations of the p-field, then we are 
left with the O(N) symmetric nonlinear sigma-model in 
the usual parametrization. 



2. The low— temperature fixed point 

Now we define the scale dimensions of the fields TTi and 
Ap as 



[n i (x)] = ~(d-2 + r,') 



[Ap(x) 



1 



[d-2 + rf) 



(B6a) 



(B6b) 



and perform a momentum-shell RG procedure.0 Here 
we have defined the scale dimension of a length L to 
be [L] = —1, and the above relations define the expo- 
nents rj and 77'. For 77, this definition coincides with 
that of the exponent usually denoted by this symbol. 
The stable, low-temperature FP of the action describes 
the ordered phase. Physically, one expects short-ranged 
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Ap-correlations, and diffusive 7r-correlations at this FP, 
which suggests the choiceS 



T] = 2 



7/ = 



(B7) 



Because of the positive scale dimension of ff, it is conve- 
nient to expand the sigma-model action, 



S 



NLoM — 



- t J dx(W(x)) 2 + ^ J dx^(x) 
+v J dx (vf(x)VTf(x)) 2 + 0(vr 4 ,V 2 7r 



(B8) 



The FP a ction is given by the first terms in Eqs. (BS) 
and (B5c), respectively. With respect to this FP, the 
magnetic field is the only relevant operator, with a scale 
dimension [h] = 2. t and r are marginal, and all other 
operators are irrelevant with respect to this FP, so the 
FP is really stable. 

The leading correction to scaling near d — 2 is given 
by the operator v, with a scale dimension [v] = — (d— 2). 
This identifies d~ = 2 as the lower critical dimension 
of the problem, v also provides the leading nonanalytic 
correction to scaling in all dimensions, since the lead- 
ing irrelevant operator related to the massive Ap-field 
has a scale dimension [c*- 2 - 1 ] — — 2, and therefore leads to 
analytic corrections to scaling. For d > d~ , the sponta- 
neously broken O(N) symmetry leads to N— 1 Goldstone 
modes, which are represented by the correlations of the 
n—1 7r-fields, 



(7r i (k)7r J -(-k))=5«|k| 



= %/k 2 



(B9) 



The FP value of the magnetization, m = (<pi (x)) , is equal 
to M, and the leading correction is given by the tt-tt 
correlation function. By virtue of the scale dimensions of 
7r and h, one finds for the leading field dependence of the 
magnetization 



m(h) = const. + h {d - 2)/2 



(BIO) 



which implies that the longitudinal susceptibility, xl = 
dm/dhr-r-> ft/ d ~ 4 )/ 2 , diverges in the h — ► limit for all 
d < 4.E3 Alternatively, the wave vector dependent zero- 
field susceptibility diverges in the homogeneous limit like 
Xl ~ |k| d ~ 4 . Notice that the leading nonanalytic correc- 
tions to scaling are given by the Goldstone modes, and 
are thus contained in the nonlinear sigma-model. 

The point of the above excercise, in the present con- 
text, is to demonstrate how much more one gets out of 
power counting after performing a symmetry analysis and 
separating the soft modes from the massive ones, as op- 
posed to doing the power counting on the action in the 
origina l <j> fo rmulation-EH The analogy to the procedure 
in Sec. Ill B] is as follows: Ap corresponds to the massive 
fields AP and AA, <j> corresponds to Q and S (the matrix 
theory cannot be expressed in terms of Q only), and tt 



corresponds to q. The external magnetic field h plays the 
role of the external frequency in the matrix theory, and 
the analogy with respect to the Goldstone modes and the 
scaling behavior of the order parameter is obvious. 



3. The critical fixed point 

Apart from the low-temperatur e FP discussed above, 
the nonlinear sigma-model, Eq. (B5b)j-also contains a 
critical FP, as was iiaticed by PolyakovJ^a and discussed 
by him and others Ejilil. This FP marks the instabil- 
ity of the low-temperature or broken symmetry phase; 
the O(N) symmetry is restored, and hence the fields n 
and a have the same scale dimension. 7? must then be 
dimensionless, which implies 



V 



2-d 



(Blla) 



in Eq. (B6a). The exponent r\ in Eq. (B6t), on the other 
hand, is the critical correlation function exponent that 
is related to the order parameter exponent /3, and the 
susceptibility exponent 7, via 



7? = 2- 7/zv = 2 -d + 2/3/v 



(Bllb) 



At criticality, the Ap fluctuations become massless, 
and the region of validity of the sigma-model shrinks 
to zero. Specifically, the Ap fluctuations can be ne- 
glected only if one works at momenta p that are larger 
then the mass of the Ap-field, m p = ^r/c^ 2 \ see Eq. 
(B5c). Power counting, and Eq. (B6b), shows that the 



scale dimension of m p is [m p ] = 1 , so upon approaching 
criticality it vanishes like m p ~ ~ \T — T c \" . The 
criterion for the validity of the sigma-model is therefore 



V <. \T - T c 



(B12) 



In this region perturbation theory works, and the critical 
properties can be explored by using the RG to sum the 
perturbative results. In this way one obtains a descrip- 
tion of the Heisenberg critical behavior in the vicinity 
of d = 2 that complements the perturbativ e trea tment 
of the 4 -formulation of the problem, Eq. (Bla), near 
d = 4. The critica l FP of the matrix nonlinear sigma- 
model, Sec. [II B 3, is again in many respects analogous 



to the critical FP of the O(N) model. 

It has been suggested that terms with more than two 
gradients, which appear as corrections to the sigma- 
model upon explicitly integrating out the massive field, 
may lead to an instability of the critical FP, even though 
thesej— terms are irrelevant by power counting at tree 
level£3. The status of this problem is currently unclear, 
see Ref. [73| for a recent discussion. 



APPENDIX C: THE INVARIANT MEASURE I[P] 

Here we derive explicitly the invariant measure I[P], 
defined by Eq. Q3.25| ). Let us write Eq. Q3.24J ) in block 
matrix form, 
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Q K> 




(CI) 



Here each block is a ANn x ANn matrix. P is block diag- 
onal by construction, and the blocks of S are subject to 
constraints due to the definition of the coset space. Now 
consider a variation SQ of Q. Since by symmetry the 
invariant measure is independent of the group parame- 
ters, it suffices to consider an infinitesimal variation. To 
first order in infinitesimal qua ntitie s, the variation of S 
is block off-diagonal (see Eq. (3.48)), and we have, 



SQ» SQ>< 
SQ<> 6Q« 

SP> SS>< P< - P> SS>< 

SS<> P> - P< SS<> 5P< 



(C2) 



By directl y d ifferentiating Q with respect to P we find 
from Eq. @, 



I[P] = det [1 ® (P<) T - P> ® 1] 



(C3) 



P > and P < are analytic continuations of anti-hermitian 
matrices, so they can be diagonalized by means of ana- 
lytic continuations of unitary ANn x ANn matrices. Let 
their eigenvalues be Xf and Xf (i = 1, . . . ANn). Putting 
the space dependence back in, we finally obtain 



I[P] 



4Nn 



n n 



(*)) 



(C4) 
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